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ABSTRACT 

We present multiple-epoch photometric monitoring in the J, H, and Kg bands of the T1.5 dwarf 
2MASS J21392676+0220226 (2M2139), revealing persistent, periodic (P = 7.721±0.005 hr) variabil- 
ity with a peak-to-pcak amplitude as high as 26% in the J-band. The light curve shape varies on a 
timescale of days, suggesting that evolving atmospheric cloud features are responsible. Using inter- 
polations between model atmospheres with differing cloud thicknesses to represent a heterogeneous 
surface, we find that the multi-wavelength variations and the near-infrared spectrum of 2M2139 can 
be reproduced by either (l)cool, thick cloud features sitting above a thinner cloud layer, or (2)warm 
regions of low condensate opacity in an otherwise cloudy atmosphere, possibly indicating the presence 
of holes or breaks in the cloud layer. We find that temperature contrasts between thick and thin cloud 
patches must be greater than 175 K and as high as 425 K. Wc also consider whether the observed 
variability could arise from an interacting binary system, but this scenario is ruled out. 2M2139 joins 
the T2.5 dwarf SIMP0136 discovered by Artigau and coworkers as the second L/T transition brown 
dwarf to display large- amplitude variability on rotational timescales, suggesting that the fragmenta- 
tion of dust clouds at the L/T transition may contribute to the abrupt decline in condensate opacity 
and J-band brightening observed to occur over this regime. 

Subject headings: brown dwarfs: general — brown dwarfs: individual(2MASS J21392676-I-0220226) 



1. INTRODUCTION 

With observed temperatures ranging from ~2200- 
300 K, field brown dwarfs (BDs) with L and T spectral 
types possess the coolest atmospheres currently available 
to direct and detailed study, and thus constitute impor- 
tant precursors to investigating directly imaged extra- 
solar giant planet atmospheres. At these cool tempera- 
tures, an understanding of atmospheric condensates — 
which affect luminosities, radii, cooling histories, and 
spectral morphologies — is essential to deriving correct 
phys ical properties for both cool BDs and g iant planets 
(e.g. iMarlev et aL|[200l iBurrows et all [20111) . 

Our current understanding of ultracool atmospheres, 
including the formation and sedimentation of conden- 
sate clouds, has developed based on comparisons of de- 
tailed atmosphere models to observations of hundreds of 
L and T dwarf s ident i fied in the solar neighb orhood (e.g.. 
Gushing et all l2008t [Stephens etldl 120091: iWitte et alj 
201 If ). Nonetheless, fundamental questions remain con- 
cerning the most basic properties of condensates, includ- 
ing their vertical and horizontal distributions, how these 
evolve as a function of effective temperature, and the 
role of secondary parameters such as gravity, metallicity, 
convection, and rotation. 

Differences in condensate clouds are responsible for a 
wide range of spectral morphologies at a given effective 
temperature, as illustrated by the transition between 
cloudy L and clear T spectral types. A detailed re- 
view of the L and T spectral sequence can be found in 
iKirkpatrickl ()2005[) and references therein. In short, the 
L spectral sequence is characterized by the formation and 



thickening of condensates with decreasing effective tem- 
perature, the L/T transition (~L7-T4 spectral types) by 
a significant decrease in condens ate opacity at a roughly 
constant eff ective temperature (jGolimowski et aP 120041 : 
[Stephens et al. 20Q9i), and the mid-to-late T spectral se- 
quence by the cooling of a cloud-free atmosphere, in 
which condensates have gravitationally settled below the 
photosphere. 

Whereas the predominantly cloudy and predomi- 
nantly clear end-points of the L/T t ransition have 
been well-modeled by various groups ([Ghabrier et al.l 
200nl: [AlIaTd et all [2?mt [fs^ [2(M lAU^^ et al.l l2003r 
Burrows et al. 1 120061 ISaumon fc Marlevll20Q8D . the tran- 
sition regime remains poorly understood. While a pro- 
gressive reddening in NIR J — Ks color is seen through- 
out most of the L spectral sequence due to an increas- 
ing photosphcric dust opacity, the L/T transition is 
evident as an abrupt blue-ward shift in J — Kg color 
by ~2 magnitudes, accompanied by a brighten i ng in 
the J-band of up to ^1 m ag (e.g. [Vrba et all \2004 . 
[Zapatero Osorio et all [20071 Faherty et al. 2011, sub- 
mitted). This dramatic blue- ward evolution and J-band 
brightening are not well-reproduced by ID atmosphere 
models inc orpo rating condensate formation and sedi- 
mentation ([Ackc rman k Marlevl[200l[:IAllard et al.l l2003l: 
[Burrows et al.|[2006c iSaumon fc Marlev[[2008D which pre- 
dict a much more gradual turn around of J — 
color with decreasing effective temperature. Vari- 
ous ideas have been put forward to explain the L/T 
transition including a sudden rain-out of condensates 
or equivalently an evolving sedimentation efficiency) 
Knapp et all [200l . or a fragmentation of the cloud 
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layer (jAckerman fc Marlevl[200lt iBurgasser et al.ll2002[ ). 
Both scenarios can reproduce main features of the tran- 
sition such as J-band brightening or the re-emergence of 
molecular FcH in early T dwarfs, which becomes depicted 
above the cloud layer due to condensate formation. 

However, neither idea has been modeled in detail as 
a physical process: both simply represent requirements 
that, if satisfied, could reproduce the observed behav- 
ior of BDs on a color-magnitude diagram. In this re- 
spect, the cloud fragmentation hypothesis is attractive 
due to its ability to make testable predictions: cloud 
structures, if possessing large enough azimuthal asymme- 
tries, should give rise to photometric variability on rota- 
tional timescales. Furthermore, in om own Solar System 
clouds on giant planets form discrete patterns coupled to 
rotationally-dominatcd atmospheric flows. For Jupiter, 
clearings in the cloud layer are evident as 5 fim hot-spots, 
and may res ult in rotationally modu lated variability as 
high as 20% (IGelino fc Marlevll2000l) . 

The combination of condensate clouds and rapid 
rotation has long-motivated searches for cloud and 
weather-related varia bility of ultracool dwarfs (e.g. 
iTinnev fc Tollevlll999f). Observations in the red optical 



jBailcr- Jones & Mimdt| |200H IGelino et all l2002t iKoenI 
[200a. ,2004, ,2005aiabl ) have targeted mainly late-M and 
early-L dwarfs due to a drop-off in optical flux for later 
spectral types. Evidence for periodic variability is found 
for ~30% of early-type objects monitored in the /-band 
with peak-to-pcak amplitudes of a few percent. How- 
ever, it is unclear whether modulation is the result of 
magnetic spots, dust clouds or a combination of both. In 
order to study variability at the L/T transition, a move 
to infrared wavelengths wherein late-L and T dwarfs are 
bright is required. 

Since the atmospheres of late-L and T dwarfs are 
increasingly neutral, they are less likely to support 
cool magnetic spots (IGelino et al.ll200^ iMohantv et al.l 
l2002f ). making the interpretation of detected variability 
in this regime less ambiguous. However, in contrast 
to optical surveys, ultracool dwarf variability at near 



and mid infrared wavel e ngths (e.g. lEnoch et al 
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Goldman et all 120081: iBailer-JonesI 120081: iClarke et alj 
2008[ ) appears to be rare (at least at sini i lar ain plitudes). 



In a study of 18 L and T dwarfs, iKoenI (|2004D found no 
significant evidence of variability in the J band above the 
~20 mmag level nor in the H or Kg bands above ~40 
mmag . A further study focussing on known /-band vari- 
ables ()Koen et al.ll2005l ) failed to detect any variability 
in the near infrared (NIR) J, //, a nd Kg bands i n any 
of these same objects. In contrast, lEnoch et al.l (|2003D 
monitored 9 objects with L2-T5 spectral types in the Kg 
band and claim large-amplitude periodic variability (0.1- 
0.4 mag) for 3 objects in their sample. However, these 
detections have not been repeated at later epochs, and 
their amplitudes are only ^1-3 times the reported pho- 
tometric precision. 

Wh ile detections at the level claimed by lEnoch et al.l 
(|2003| ) have not been reported since, recent results from 
high-precision NIR monitoring provide some evidence for 
lower-amplitude variability in the N IR. In a surv e y of 8 
late-L and T dwarfs in the J band, IClarke et ahl (|2008D 
claim periodic variability for 2 of 7 late-L and T dwarfs 
monitored with amplitudes of 15 and 8 mmag and pe- 



riods of 1.4 hr and 2 hr r espectively. Most r ecently, 
in a single targeted study, lArtigau et al.l (I2009D found 
the T2.5 dwarf, SIMP0136 (|Artigau et al]l2006D . to be 
variable with a peak-to-peak amplitude of ^--^50 mmag in 
J and a period of 2.3 hr (a 10-cr detection). The high 
significance, repeatability and multi-band detections for 
SIMP0136 represents a breakthrough in comparison to 
previous work. Notably, SIMP0136 falls directly within 
the L/T transition regime. 

Here we present continuous photometric monitoring 
in the J, //, and Kg bands of the Tl 5 dwarf 2MASS 
J21392676+0220226 (|R,eid et al.l [20081 2M2139 here- 
after) using the Wide Field Infrared Camera on the 
DuPont 2.5-m telescope at Las Campanas. 2M2139 is the 
most variable BD to date, and was discovered in our large 
J-band monitoring program for weather-related variabil- 
ity of cool BDs (Radigan et al., in preparation). In S}2] 
we describe the observations and data reduction, and in 
fJSlwe present and analyze the reduced light curves. In 
f|4lwe present an overview of 2M2139 's physical proper- 
ties based on existing archival data. In ^we attempt to 
model the observed multi-band variability using linear 
combinations of 1-dimensional cloudy and clear atmo- 
sphere models. In fj6]we discuss our results and possible 
origins of the observed variability. Finally, in SJ7] we give 
a brief summary of our conclusions and suggestions on 
how to proceed. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. LCO Observations 

Observations of 2M2139 were made using the W ide 
Field Infrared Camera (WIRC; iPersson et aLll2002D on 
the Du Pont 2.5 m telescope at Las Campanas, as part of 
a large survey for J-band variability of cool BDs (Radi- 
gan et al., in prep). The camera consists of 4 HAWAII-I 
arrays, each with a 3.2' field of view and a pixel scale of 
0.2". The camera is intended as a wide-field survey cam- 
era, with 3' gaps between detectors. We have not used 
it as such, choosing to position our target consistently 
on the south-west array, which we determined to be the 
least noisy of the four chips. New J-band filters closely 
matching the Mauna Kea Observatory (MKO) system 
(jTokunaga fc Vaccal [20051) were purchased and installed 
for our survey in order to minimize the effects of differen- 
tial atmospheric extinction, by cutting off time-variable 
telluric water absorption bands red- ward of 1.35 /im (e.g, 
lArtigaul [20061 ). 

2M2139 was first monitored in the J band on 01 Aug 
2009 for 2.5 hours. An analysis of its light curve re- 
vealed a 90 mmag increase in brightness over this time 
span. This prompted us to follow-up with longer pho- 
tometric sequences in J. H, and Kg bands on the four 
consecutive nights of 21-24 Sep 2009. In order to obtain 
near-simultaneous photometry in multiple bandpasses, 
additional photometric sequences alternating between J, 
//, and Kg bands every ~20 min were obtained on 26 
Sep 2009, 30 Sep 2009, and 01 Oct 2009. A log of all 
observations is provided in Table [TJ 

All observations except those of 21 Sep 2009 were made 
using a random dither pattern wherein the telescope was 
offset by at least 3"(15 pixels) after each exposure, and 
all pointings were contained in a 15" (75 pixels) square 
box. Individual exposure times in J were 40 s (09 Aug 
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Table 1 

Observing Log 





Instr.'^ 


Filter 


UT Start'' 


Length (hr) 






Notes 




01 Aug 2009 


WIRC 


J 


00hl9m 


3.1 


40 


140 


Clear at start with light clouds after 1.5 hr, 0.7' 


seeing 


21 Sep 2009 


WIRC 


J 


00hl4m 


5.12 


60 


255 


Clear with occasional passing light cloud, 0.7" 


seeing 


22 Sep 2009 


WIRC 




23h53m 


5.48 


20 


412 


Clear with occasional passing light cloud, 0.8 


seeing 


23 Sep 2009 


WIRC 


J 


23h29m 


5.99 


60 


260 


Clear with occasional passing light cloud, 0.7 


seeing 


24 Sep 2009 


WIRC 


H 


23h25m 


6.01 


20 


498 


Clear with occasional passing light cloud, 0.7 


seeing 


26 Sep 2009 


WIRC 


J,H,K, 


23h45m 


5.58 


45,20,20 


91,130,132 


Clear, 1.5" seeing 




30 Sep 2009 


WIRC 


J,H,K^ 


23h40m 


4.81 


45,20,20 


82,110,114 


Clear at start becoming cloudy after first 4.1 hr, 1.2"socing 


01 Oct 2009 


WIRC 


J,H,Ks 


23h46m 


5.11 


45,20,20 


77,104,110 


Clear, 1.1 seeing 




08 Nov 2009 


CPAPIR 


J 


22h37m 


4.37 


21.6 


387 


Clear, 1.5-2" seeing 





Dates correspond to the local day at the beginning of the night. 

WIRC observations wore conducted on the LCO 2.5-m telescope, and CPAPIR on the OMM 1.6-ni telescope 
UT Start time may be +1 day ahead of the local date. 



2009), 45 s (26 Sep 2009, 3 Sep 2009, 01 Oct 2009), or 
60s (23 Sep 2009). Exposures were read out using cor- 
related double sampling. Along with telescope offsets 
this resulted in effective cadences of 62 s, 68 s, and 82 s 
respectively for the listed exposure times. For the H- 
and Xs-bands we used individual exposure times of 20 s 
throughout, resulting in an effective cadence of 42 s. 

The J-band observations of 21 Sep 2009 employed 
a staring strategy, in which the target centroid was 
kept fixed on the same pixel throughout the sequence. 
This was accomplished using an IDL routine to stream 
the incoming science images onto a standard laptop 
and compute real-time guiding corrections. An alert 
was sounded, and manual closed-loop corrections to the 
guide-camera reference position were made each time the 
target strayed by more than 0.5 pixels from its initial po- 
sition (approximately once every 5-15 minutes). Without 
telescope offsets the efficiency of staring observations is 
significantly increased. Individual exposures of 60 s were 
used, resulting in a cadence of 67 s. Nine-point dither 
sequences for the purpose of rough sky subtraction and 
centroiding were made at the beginning, middle and end 
of the 5.12 hr observation. 

For each sequence, dome-flats (lamp on and off) and 
dark frames corresponding to each exposure time were 
taken cither on the afternoon preceding, or the morning 
following each observation. 

2.2. Reduction and Processing of WIRC data 

All raw images were corrected for non-linearity using 
a detector response curve measured on 27 Jul 2009. Cal- 
ibration images were median-combined to create high 
signal-to-noise dark and flat-field frames. The dark- 
current contribution was then subtracted from all other 
images. For the dithered sequences, a running sky frame 
was computed by median combining the 11(7) J{HK)- 
band images that were (i)taken closest in time to, and 
(ii)were spatially offset by at least 6" (30 pixels) from 
the image being reduced. After an initial first-pass reduc- 
tion, stars were identified from the stacked field, and then 
masked for the second pass so as not to bias the median- 
combined sky frames. After sky subtraction, pixel-to- 
pixel variations in quantum efficiency were removed by 
dividing the resultant image by the fiat-field. Bad pixels 
identified from the flat-field or pixels having more than 
35000 counts (>3% non-linear) were flagged. Except for 
sky subtraction, reduction of the non-dithered staring 
sequence is almost identical to the procedure described 



above. In this case the sky frames were constructed by 
median combining the 9-point dither sequences taken 
before and after staring, linearly interpolating between 
them in time, and scaling to the 3cr-clipped median of 
each science image. Since immediate background lev- 
els are also subtracted using aperture photometry, the 
primary motivation of sky subtraction for the staring se- 
quence is to achieve better centroiding of the target and 
reference PSFs. For all reductions and analyses we con- 
sidered only the south-west array of the WIRC camera. 

2.2.1. Relative Differential Ptiotometry 

For each monitoring sequence, aperture photometry 
was performed on 2M2139 and a set of reference stars in 
2M2139's the field of view, using a circular aperture of 1.5 
times the median full width at half maximum (FWHM) 
of all stars in each image. A large aperture was chosen 
to decrease the systematic effects of slightly elongated 
PSFs on our photometry which occasionally occur due 
to a gradual degradation of the telescope focus during 
long sequences. Residual sky levels in the vicinity of each 
star were measured inside an annulus centered on each 
source of inner radius 3.7 times the FWHM and a width 
of 11 pixels. Flux measurements in which a fiagged pixel 
fell inside the aperture were set to an error value and 
disregarded in the light curve analysis. 

The raw light curves display fiuctuations in brightness 
due to changing atmospheric transparency, airmass, and 
residual instrumental effects throughout the night. To 
a very good approximation these changes are common 
to all stars, and can be removed. First, the raw light 
curves of all stars were converted from absolute to rela- 
tive fluxes via division by their median brightness. Next, 
for each reference star a calibration curve was created by 
median combining the light curves of all other reference 
stars (excluding that of the target and star in question). 
The raw light curve of each reference star was then di- 
vided by the corresponding calibration curve to obtain a 
corrected light curve. The standard deviations of the cor- 
rected light curves for each reference star (cr) were then 
measured. This process was repeated multiple times, us- 
ing an iterative approach where only reference stars with 
a < 15mmag in J and <25mmag in H and Ks were 
keptQ This process was terminated (usually after the 
second iteration) when the number of good calibration 

^ For the observations of 30 Sep 2009 these criteria were relaxed 
to fT < 20 mmag in J and <30 mmag in Kg , and for the observations 
of 01 Oct 2009 they were relaxed to a OOmmag in Kb- 
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Table 2 

Reference stars used for differential pliotometry 



2MASS ID 


Letter ID" 


J 


J-K, 


Epoelis used ^ 


J21392216+0220185 


A 


15.51 


0.65 


2,4,9 


,121392311+0222009 


B 


14.46 


0.64 


2-9"^ 


J21392542+0222102 


C 


12.98 


0.48 


1-3,5-9 


J21392392+0222383 


D 


15.22 


0.93 


2-5,6-7(J, H only), 8-9 


J21392228+0223082 


E 


14.26 


0.64 


2,4-7 


J21392465+0223140 


F 


14.63 


0.92 


2-9 


J21393173+0222126 


G 


11.48 


0.37 


3,5-9 


J21393318+0222356 


H 


15.40 


0.52 


1 


J21393533+0220584 


I 


13.50 


0.36 


1"= 


J21393502+0220466 


J 


15.33 


0.59 


1 


J21393596+0220488 


K 


15.65 


0.59 


1 



" stars are labelled by letter in figure [Tl 

^ Epochs from tableware numbered chronologically from 1-9. 
Comparison star 



stars no longer changed. 

A calibration curve for 2M2139 was then computed us- 
ing this subset of high signal-to-noise (S/N) references, 
loss a comparison star of similar brightness. Figure [T] 
shows a stacked and mosaicked image of the field with 
reference stars labelled by letter. Details (identifiers and 
magnitudes) of the reference and comparison stars used 
to calibrate each set of observations are provided in ta- 
ble [5] We note that on some nights otherwise high S/N 
reference stars were excluded due to consistently falling 
on bad pixels, saturation, or contamination from a faint 
nearby source in poor seeing conditions (e.g. star E). 
While it would be optimal to use a common set of refer- 
ence stars throughout, we find that the shape and am- 
plitude of our final light curves are independent of our 
choice of reference stars. For most of our observations 
star B from figure [1] was chosen as the comparison, with 
the exception of 01 Aug 2009 when it didn't fall in the 
field of view and we used star / instead. 

Detrended light curves for 2M2139 and the comparison 
star were obtained by dividing their raw fiuxes by the 
final calibration curve. An example of raw and detrended 
light curves for the reference stars and target, as well 
as light curve standard deviations as a function of star 
brightness for the 23 Sep 2009 epoch are shown in figure 
H Detrended J, H and light curves for 2M2139 and 
the comparison star for the consecutive nights of 21-24 
Sep 2009 are displayed in figure |31 

In addition, light curve properties for all epochs in- 
cluding amplitudes, and standard deviations for the tar- 
get (at) and comparison star (cc) are provided in table[3l 
Since the light curves for 2M2139 show large trends, at is 
obtained by taking the standard deviation of 2M2139 's 
light curve subtracted by a shifted version (by one ele- 
ment) of itself, and then divided by \/2. For the epochs 
where we cycled through the J, H, and Kg filters, large 
time gaps exist in the sequence for any given filter, ne- 
cessitating a different method for obtaining at- Thus 
for these sequences we took at to be the average of the 
standard deviations obtained for individual ^20 min seg- 
ments in a given filter. 

2.3. OMM Observations 

Due to the evolving nature of the light curve, the tar- 
get was observed again on the night of 08 Nov 2009 using 
the Observatoire du Mont Mcgantic 1.6-m telescope and 
the Camera PAnoramic Proche Infra-Rouge (CPAPIR, 
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Figure 1. Stacked and mosaicked WIRC J band image of the 
field surrounding 2M2139 . The target and reference stars used for 
differential photometry are labeled, with details provided in table 

m 

lArtigau et aT1[200l . The camera has a 30' field of view 
with a pixel scale of 0.89". As with our WIRC obser- 
vations, an MKO J filter was employed. Observations 
commenced at an airmass of 1.39 and concluded at an 
airmass of 2.35. Random dithers were made between 
each exposure, always keeping the target within 1' of the 
central position. A single co-addition with an exposure 
time of 21.6 s was used for all science frames, resulting in 
a median cadence of 38 s. 

2.4. Reduction of CPAPIR data 

The raw CPAPIR images were processed using the 
CPAPIR pipeline. A running median sky image were 
constructed by taking the median of the 11 frames 
centered on the frame of interest after masking bright 
field star s identified in the Tw o Micron All Sky Survey 
(2MASS: ISkrutskie et ai1l2006[) . After sky subtraction, 
images are divided by a flat fleld image constructed from 
on/off dome images. The Poisson noise associated with 
the fiat field is ^0.2% per pixel and a negligible contrib- 
utor to the overall error budget. Aperture photometry 
was performed in a similar manner to the WIRC images, 
but using an aperture size of 1 FWHM, and an annulus 
of inner radius 3.5 FWHM, and 11 pixel width. 

2.5. Photometric Calibration 

In addition to the relative photometry, all light curves 
were flux-calibrated against the 2MASS catalog using ref- 
erences identified in figure [T] While the WIRC H and Kg 
fihers are similar to those from 2MASS, the MKO J fil- 
ter installed for our survey is significantly narrower than 
the 2MASS J filter. Because we do not have enough ref- 
erence stars to derive robust color-dependent corrections 
between our MKO J filter and 2MASS J filter, we first 
converted 2MASS J magnitudes of the reference stars 
to MKO magnitudes accord ing to the color-dep endent 
transformations provided bv lLeggett et al.l ()2006[ ). 

The 2MASS catalog magnitudes, corrected for our fil- 
ter system, were then converted to counts (multiplied by 
some arbitrary factor) and a linear fit between the cat- 
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Figure 2. Left: An example of raw (top) and detrended (bottom) liglitcurves for 2M2139 (grey lines) and reference stars (colored lines) 
for a photometric sequence taken on 23 Sep 2009. Right: Standard deviations of the detrended light curves for all stars on the chip as a 
function of J magnitude. 2M2139 is shown as an open star symbol. A solid line shows the photon noise limit. At J magnitudes brighter 
than ~16.5 observations are close to being photon-limited. 



alog and measured counts was performed. The quality 
of the photometric calibration was determined from the 
slope error for each linear fit, and found to be accurate 
among reference stars at the 1-2% level. 

In figures we have opted to plot light curves in units 
of relative flux as opposed to magnitudes. The refer- 
ence magnitudes corresponding to relative fluxes of 1 are 
given by Jmto=14.75, H=U.n, and K,=U.59. The 
chosen reference magnitudes correspond to the median 
magnitudes measured in each bandpass for the near- 
simultaneous JHKg sequence of 26 Sep 2009, where the 
light curves were observed to have similar amplitudes in 
all bandpasses. 

In addition to placing all light curves on a common 
flux scale, the photometric calibration allows us to com- 
pare the magnitudes and colors derived here with archival 
data from 2MASS and the SpeX Prism Librarjfl. In the 
SpeX Prism Library there are two low resolution NIR 
spe ctra (R~120) for 2M21 39 , obtained in 2003 and 2004 
by iBurgasser etall (|2006[ ). The SpeX Prism Library, 
while not able to provide absolute photometry, allows us 
to derive colors for 2M2139 at additional epochs. When 
synthetic 2MASS colors for M, L and T dwarfs from the 
SpeX Prism Library are compared to their 2MASS val- 
ues, the match is surprisingly good, with a standard devi- 
ation similar to the reported photometric errors from the 
2MASS catalog. For L and T dwarfs no significant sys- 
tematic trends exist as a function spectral type, whereas 
we find the synthetic J — Kg colors of M dwarfs to be 
on average slightly bluer than values derived in 2MASS. 
We present a comparison between 2MASS and synthetic 
SpeX colors in Appendix [X] 

Finally, for a proper comparison between WIRC, 
2MASS, and SpeX (synthetic 2MASS) photometry, we 
have converted WIRC magnitudes derived for 2M2139 to 
the 2MASS system. Correction terms between the two 
filter systems were computed directly from 2M2139's 
NIR spectrum using 2MASS and WIRC filter plus sys- 
tem transmission curves and a Kurucz mo del Vega spec- 
trum as a zero-magnitude fiux reference (jKuruczl 119791 
|1993[) . Recall that for reference stars corrections between 

^ http:/ /pono. ucsd.edu/~adam/browndwarfs/spexprism/ 



J filters were accounted for in advance by our conver- 
sion of 2MASS to MKO J magnitudes, while H and Ks 
band corrections are generally negligible. Magnitudes 
and colors measured for 2M2139 from WIRC, 2MASS 
and SpeX epochs, as well as correction terms specific to 
2M2139 used to convert between WIRC and 2MASS fil- 
ter systems, are provided in table S] 

3. RESULTS AND ANALYSIS 

3.1. The Reduced Light Curves 

Full J; i?, and Kg light curves from the adjacent nights 
of 21-24 Sep 2009 are shown in figure|3l These epochs en- 
compass the largest variability observed for 2M2139 with 
peak-to-peak amplitudes of 26% in J, 17% in H and 16% 
in Ks- Peak-to-peak amplitudes are measured directly 
as the maximum change in brightness, divided by the 
mid-brightness. We note that the 77-band variation may 
be slightly larger since its global maximum is not entirely 
captured. 

The near-simultaneous J. H , and Kg light curves taken 
on 26 Sep 2009, 30 Sep 2009, and 01 Oct 2009 are 
shown in figure 21 Unbinned data are shown in the 
top panel, while in the middle panel the data have been 
binned to one data point per filter change (6 exposures 
per epoch in J, and 10 in H and Kg^ less exposures 
discarded due to bad pixels). Error estimates for the 
unbinned data points are given by the standard devia- 
tion of all measurements in a continuous segment at a 
given filter position, while the binned data points were 
assigned uncertainties of 1/a//V times lower, where TV 
is the number of binned exposures at a given filter po- 
sition. Relative amplitudes of variability in the differ- 
ent bands were measured by assuming all light curves to 
have the same functional form, but with differing ampli- 
tudes. For this purpose we used the J-band light curves 
as templates, linearly interpolating between binned data 
points. Scaled versions of the J-band templates were 
then fit to the H and Kg light curves using a simple 
weighted linear regression. The best-fit solutions (which 
we will refer to as the H and Kg band templates) are 
overplotted on the data in figure|4l We find amplitude ra- 
tios of Ai//^j={0.91±0.07, 0.84±0.08, 0.91±0.15} and 
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Figure 3. Reduced light curves for 2M2139 obtained from 21-24 Sep 2009 (top) and a comparison star of similar brightness (bottom). 
Vertical lines separate observations on different nights. The given dates correspond to the beginning of each observing night. The J-band 
data points show unbinned 60s exposures, while the H and Ks data points show 20s exposures binned by a factor of 3. 



AKjAj={Qm ± 0.08, 0.59 ± 0.07, 0.45 ± 0.11} for the 
earliest to latest epoch respectively. Here, Ai represents 
a peak-to-peak amplitude of photometric variability in 
a bandpass i given by the absolute change in brightness 
divided by the mid-brightness. Uncertainties for the am- 
plitude ratios were determined by repeating the above 
procedure on a set of 1000 simulated light curves. The 
simulated light curves were obtained from sampling the 
J, -fT, and Kg band templates at the times of our binned 
data points in figure 21 A random noise component was 
added to each point, drawn from a gaussian distribu- 
tion with a standard deviation equal to the la uncer- 
tainty corresponding to that data point. In practice we 
found that the fits between the H and Kg light curves 
and scaled J-band templates produced reduced val- 
ues from 1.4-3, and we therefore scaled the width of the 
random error component to achieve a reduced ^ 1 
in comparison to the best-fit templates. Measured light 
curve amplitudes and ratios for the different epochs are 
tabulated in table [31 

The results from the near-simultaneous sequences are 
somewhat surprising. The light curves from the consecu- 
tive nights of Sep 21-24, where observations were made in 
a single but different band each night (figure [S]), would 
have led us to conclude that variations in the H and 
Ks bands are only 50-60% of those in J. However, am- 
plitude ratios measured from the first near-simultaneous 
sequence on 26 Sep 2009 appear close to unity. Upon 
repeating the near-simultaneous measurements at two 
subsequent epochs we continued to find Ah/Aj > 0.84, 
while we found Ak^/Aj as low as 0.45 in the latter epoch. 
Our data may hint at amplitude ratios that are them- 
selves variable, as amplitude ratios from the first and 
last epoch with simultaneous JHKg data are incompat- 
ible at the 2(7 level. In addition, none of the peak to 
peak amplitudes measured for the simultaneous JHKg 
light curves (irrespective of bandpass) are as large as the 
lowest amplitudes measured from the single-band con- 
secutive night sequences of 21-24 Sep 2009. Therefore, 
there is evidence that both amplitude ratios, as well as 
the overall amplitude may vary with either epoch, light 
curve phase, or a combination of both. 



3.2. Rotation and Longer Timescale Trends 

While not strictly uniform from epoch to epoch, the 
variations exhibited by 2M2139 show a clear periodic- 
ity. Since the J-band light curves do not seem to have 
evolved significantly in shape from 21 Sep 2009 to 23 
Sep 2009 we determined an approximate period by fit- 
ting a periodic function to these light curves. Since vari- 
ations are not perfectly sinusoida l, we chose a round, 
non-overlapping, two-spot model (jPorren ' 1987) as our 
periodic function, and p erformed fitting using a Markov 
Chain Monte Carlo (e.g. iLewis &: Bridld [ 2002D technique 
with a Metropolis H astings algorithm (|Metropolis et all 
119531 lHastii"Igill970[ ) (figure [Q). Using this method we 
determined a period of 7.721 ± 0.005 hr, with the quoted 
value and errors derived from the maximum likelihood 
and 67% credible regions of the posterior distribution of 
periods. The uncertainty is likely underestimated due to 
the implicit assumption that the light curve shape has 
not evolved between cycles used for fitting. 

Two spots are required to fit the asymmetry in the 
light curve, which is hinted at from a secondary "bump" 
present at the beginning of the J-band light curve from 
21 Sep 2009, and which also becomes apparent when ob- 
servations are phased to the best-fit period. 

In addition to fitting for the period, we also com- 
puted a L onib-Scarglc Pcriodogram (IScargld [l982l: 
iHorne fc Baliunas 1986) of the J-band data spanning 21 
Sep 2009 to 01 Oct 2009 (figure E]). In order not to favor 
epochs with more densely sampled data, all light curves 
were linearly resampled at 10-min spacing within each 
observing window. Despite falling close to the 8 hr sub- 
peak of the window function (which peaks at 24 hr), the 
strongest peak at 7.73 hr matches the 7.721±0.005 hr 
recovered from our spot modeling. 

We find that the uncertainty in the best-fit rotation 
period of 7.721 ± 0.005 hr gives us the freedom to phase 
the Aug 2009 light curves with those from Sep 2009, but 
there is no single period that can accommodate all epochs 
spanning Aug 2009 - Nov 2009. As an example of this, 
all light curves (presented on a common fiux scale) are 
shown in figure[71 phased to an over-constrained period of 
7.723 hr. In fine-tuning the period to precisely 7.723 hr. 
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Figure 4. Top row panels: Unbinncd interleaved J (black), H (blue), and Ks{r:cd) photometric sequences from 26, 30 Sep 2009 and 01 
Oct 2009 epochs. Unity corresponds to Jmko^ H, and Ks magnitudes of 14.75, 14.11, and 13.59 respectively. Middle row panels: Binned 
(see text for details) J, H, and Ks light curves, with same color scheme as above. Linear interpolations of the binned J band points are 
shown as black linos. Scaled, best-fitting versions of these J-band templates to the H and Ks light curves are shown as red and blue lines 
respectively. Bottom row panels: Changes in J — H and J — Ks colors, based on the scaled templates. Bottom panel: Amplitude ratios 
and uncertainties derived from Monte Carlo simulations for each of the three epochs. 



we have arbitrarily chosen to phase the earlier (Aug 2009) 
observations with those from Sep 2009 in order to illus- 
trate the resultant mismatch in phases for the later (Nov 
2009) epoch. However, given the lack of a common phas- 
ing across all epochs, there is no reason to expect the 01 
Aug 2009 observations to be in phase with the middle 
epochs. Thus, while we have used a period of precisely 
P = 7.723 hr in our figures, it should be understood 
that this is illustrative, and that the actual periodicity 
cannot be constrained beyond 7.721 ± 0.005 hr from our 
observations. 

While the maxima and minima of the light curves 
spanning 21 September 2009 to October 01 2009 appear 
roughly synchronized with a 7.723 hr period, there are 



large differences in amplitude, and more subtle differ- 
ences in light curve shape from epoch to epoch. 

There is a marked decrease in the observed light curve 
amplitudes between the 23 Sep 2009 and 26 Sep 2009 
epochs. While this may indicate short timescale evo- 
lution of cloud features, it is also consistent with a 
more stable, double-peaked light curve with a period of 
2x7.72=15.44 hr such that the higher-amplitude varia- 
tions observed on 21 Sep 2009 and 23 Sep 2009 sam- 
ple a different phase of the light curve than the lower- 
amplitude variations observed at subsequent epochs. 
This possibility is illustrated in figure [5] where light 
curves spanning 21 Sep 2009 to 01 Oct 2009 are shown 
phased to both 7.723 hr and 15.446 hr (once again, over- 
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Figure 5. Top: Best-fitting round, two-spot model (red curve, 
P=7.721 hr) to the J-band light curves from 21 and 23 Sep 2009 
(data points). 

constrained) periods. For the case where P = 15.446 hr. 
ovcraU changes in hght curve shape and amphtude, al- 
though still apparent, are less pronounced. A period of 
15.44 hr would make 2M2139 a somewhat slow rotator 
in comparison to other ultracool dwarfs with v sin i mea- 
surements, which have inferred p eriods ranging from ^2- 
12 hr Ce.g. lReiners fc Basri|[2008l ). 

There appears to be some phase coherence on a 
timescale of days to weeks, suggesting that we are observ- 
ing cloud features that persist at least this long. Over 
months the phase coherence is lost, which may indicate 
evolution of the cloud coverage or the dissipation and for- 
mation of new features. An intriguing possibility is that 
the long-timescale evolution of the light curve is caused 
by the differential rotation of a storm syste m with respect 
to clo ud features at different latitudes fe.g. lArtigau et all 
12009( 1. The light minimum of the last epoch (8 Nov 2009) 
is approximately 3 hr ahead of the continuous J- band 
sequence from 23 Sep 2009, which could imply a differ- 
ential rotation of ^3 degrees per day or wind speeds of 
~45 m during that time period (approximately half 
of that on Jupiter). We note that this estimate relies 
on our over-constrained period of 7.723 hr which was ob- 
tained by assuming a constant phase for the Aug through 
Oct epochs, which may not be realistic. Long term mon- 
itoring of 2M2139 should be able to verify the existence 
of a persistent, differentially rotating feature. 

3.3. Possible long-term variability 

From comparing our WIRC photometry to that from 
2MASS (both on the 2MASS photometric system, see 
table m, we find 2M2139 to be significantly bluer (J - 
/v 3 =1.26- 1.34) at the time of the WIRC observations 
than reported by 2MASS (J - A',.=1.67±0.07). We can 
estimate 2M2139's maximum change in J — Ks color 
within the time spanned by our WIRC observations by 
making the conservative assumption of Ak^/Aj=0A5 
(see figure 3] and table [3]). This amplitude ratio, the 
lowest in our data, can be used to infer a maximum 
A( J — Ks) of 0.15 mag for the full range of J band vari- 
ability observed (A,/=0.26). Therefore, there is no way to 
reconcile the 2MASS catalog color of J - Ks=l.67±0.07 
with that observed by WIRC, unless large systematic er- 
rors exist in our photometric calibration. Fortunately, 
NIR spectra for 2M2139 from the SpeX Prism Library 
provides two additional epochs which wc can compare to 
WIRC and 2MASS photometry. Synthetic 2MASS colors 




0.5 1.0 1.5 

Angular Frequency [rad hr '^] 

Figure 6. Lomb-Scargle Periodigram (black line) of light curves 
spanning 21 Sep 2009 to 01 Oct 2009. The power corresponding 
to a false alarm probability of 0.997 is ~4, falling well below any 
peaks in the periodigram. The scaled window function (red line) 
shows a main peak at 24 hr corresponding to the spacing of our 
nightly observations, with additional peaks at 12 hr, 8 hr, etc. The 
strongest data peak at 7.73 hr matches the 7.721±0.005 hr value 
recovered from fitting a periodic function to the light curve (see 
figure O . 

of J-Ks=l.32±0.083 and 1.36±0.080 were found for the 
2003 and 2004 SpeX epochs respectively, in agreement 
with the WIRC photometry. It is interesting to note 
that compared to 198 other L and T dwarfs for which we 
have also measured synthetic SpeX photometry (see Ap- 
pendix 2M2139's large offset between 2MASS and 
SpeX colors makes it a 4(T outlier. Thus it may have 
been possible to identify this target's variability well in 
advance of our observations. It is also notable that there 
is no apparent increase in scatter between 2MASS and 
SpeX photometry across the L/T transition, with a few 
other less extreme outliers occurring at earlier spectral 
types. These additional outliers should be followed up, 
although it appears that variability on the scale observed 
for 2M2139 is quite rare. 

Thus, if the 2MASS epoch can be trusted, there is ev- 
idence that the entire range of 2M2139's variability is 
larger than captured by our observations. We note that 
there is no reason to suspect an error in the 2MASS col- 
ors, as catalog magnitudes derived from profile fitting 
agree well with those from aperture photometry, and 
the 2MASS J, H, and Kg observations are simultane- 
ous. Furthermore, according to our WIRC photometry 
the noted blucward shift in J — Ks since the 2MASS 
epoch is almost entirely attributable to a brightening in 
the J-band, accompanied by a plateau or possible dim- 
ming in the Ks band. If real, this long-term trend of 
anti-correlated behavior in J and Ks is qualitatively dif- 
ferent than the highly correlated variability we observed 
over short timescales. 

4. PHYSICAL PROPERTIES 

The physical properties of 2M2139 derived here and 
elsewhere are presented in table |4l 2M2139 is a cool 
field BD with an optical spectral type of TO deter- 
mined by iReid et alT (I2008D, a nd a NIR spectral type 



of T1.5 from |Burgasscr et all (I2006D The NIR spec- 
tra (R'--^120) from iBurgasser et all (l2006l) were taken 
with the Medium Resolution Near Infrared Spectro- 
graph (SpcX) at the NASA's InfraRed Telescope Fal- 
cility (IRTF), and are available online from the SpeX 
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Table 3 

Light Curve Properties 



Epoch^ 


Filter 


A" 


Ah/Aj 


Ak, /Aj 






01 Aug 2008 


J 


0.080 






0.0072 


0.0056 


21 Sep 2009 


J 


0.260 






0.0067 


0.0051 


22 Sep 2009 


Ks 


0.156 






0.0145 


0.0156 


23 Sep 2009 


J 


0.262 






0.0079 


0.0070 


24 Sep 2009 


H 


0.166 






0.0124 


0.0119 


26 Sep 2009 


J {H,K^) 


0.134 (0.123, 0.111) 


0.91±0.07 


0.83±0.08 


0.0092 (0.0159,0.0186) 


0.0088 (0.0136,0.0236) 


30 Sep 2009 


J iH,Ks) 


0.109 (0.092, 0.064) 


0.84±0.08 


0.59±0.07 


0.0107 (0.0154,0.0201) 


0.0107 (0.0132,0.0227) 


01 Oct 2009 


J (H,Ks) 


0.084 (0.077, 0.038) 


0.91±0.15 


0.45±0.11 


0.0088 (0.0142,0.0175) 


0.0076 (0.0116,0.0188) 


08 Nov 2009 


J 


0.152 






0.0391 


0.0232 



^ As in table ^ dates correspond to the local day at the beginning of the night. 

^ Peak-to-peak amplitudes arc measured as the absolute change in flux, divided by the mid-flux and correspond to the filtcr(s) indicated in the 
second column. 



Prism Library. There is no parallax data available for 
this target at this time, and given it's unique vari- 
ability (and hence potentially unique atmospheric char- 
acteristics) it is unclear to what degree standard re- 
lationships between spectral type, absolute magnitude 
and temperature may apply. For field BDs (~3 Gyr) 
with measured parallaxes and estimated bolometric lu- 
minosities it has been shown that effective temperatures 
plateau or slightly decrease from ~1 400-1200 K across 
the L/T transition ( ^L7-T 4 SpTs; iGolimowski et alj 
120041: istephens et"all I2OO90 . The scatter about this 
temperature is repor ted as ^100 K . Using the rela- 
tionship provided by [Stephens et al.l ()200 9) we derive 
Te//=1270±100 K for a T1.5 spectral type. However, the 
temperature of the L/T transition also shows some de- 
pendence on gravity fe.g. iMetch ev & Hi llenbrandl 120061 : 
iLeggett et al.ll2008HStephens et al.,.2009l and this tem- 
perature could be lower if 2M2139 is moderately young. 

If we assume 2M2139 to have an age/mass typical of 
field T dwarfs, an absolute magnitude and distance can 
be determined from empirical relationships. Using the 
sp ectral type versus MK O K magnitude relation given 
bv iMarocco et al.l (|2010[ ). excluding known binaries and 
using a NIR spectral type of T1.5, we find an absolute 
magnitude of Mk — 13.35 ± 0.25 mag. The error bar 
refiects the coeffici ent uncer tainties of the polynomial fit 
provided by Maro cco et al.l (j2010), but neglects intrinsic 
scatter of the sample about the relati on. This intrinsic 
scatter is estimated bv lLiu et al.l ()2006[ ) to be ~0.39 mag, 
and we therefore adopt this larger uncertainty hereafter. 
Assuming 2M2139 is a single object, this corresponds 
to a distance of ILG^i gPC, where we have used the 
2MASS Kg magnitude of 13.58 mag, first converted to 
an MKO K magnitu d e usin g the relationship provided by 
Istephens fc LeggettI (|2004l ). to obtain the distance mod- 
ulus. Our distance esti mate is notably closer than that 
of 18.8 pc cataloged b v iLooper et al.l (|2007t ) using the 
MKO- J relationship of lLiu et al.l ()2006t ). 

4.1. Binarity? 

iBurgasser et al.l (|2010f ) have suggested that 2M2139 is 
a binary candidate since a composite spectral template 
consisting of L8.5±0.7 and T4.5±1.5 components pro- 
vide a significantly better fit to its NIR spectrum than 
any single template. For example, pronounced CH4 ab- 
sorption is present at 1.2 /xm, but is weak or absent in H 
and Ks ■ One way to reproduce this mismatch is with the 
addition of a late T companion with highly suppressed 
H and Kg fluxes; this would enhance the height of the J- 



Table 4 

Target Properties 



Quantity 


Value or Range 


Data Rcf. 


Source 


Identifier 


J21392676-I-0220226 


2MASS° 


2 


a (J2000) 


21 '"39™ 26". 76 


2MASS'' 


2 


<5 (J2000) 


H-02'' 20'22".6 


2MASS'' 


2 


Optical SpT 


TO 




3 


NIR SpT 


T1.5 




4 


Period 


7.721±0.005 hr 




1 


Te// (3 Gyr) 


1270 ± 100 K ° 




1 


Mk.mko (3 Gyr) 


13.35 ± 0.25 ' 




1 


d 






1 


Photometry 


,7 


15.26±0.049 






H 


14.16±0.053 


2MASS'' 


2 




13.58±0.045 






J - 


1.67±0.066 






J - H 


1.10±0.072 






J - K, 


1.32±0.07 






J - H 


0.89±0.07 






^J{2M-MKO) 


0.183 


SpeX'= 


1 




-0.065 








-0.016 






'HiKO 


14.66-14.96'= 






J 


14.84-15.14"= 






H 


14.00-14.18'= 


WIRC 


1 


Ks 


13.54-13.71"= 






J - Ks 


1.26-1.34'' 






J - H 


0.78-0.83"* 







No te. — References r efer to (l)This paper. (2lSkrutskie et all 
120061) . (3 1Reid et al.l 120081) . (4 1Burgasser et al.l II2006I ) 
" 2MASS Point Source Catalog, the epoch is JD 2451741.8470 (16 Jul 
2000) 

^ Synthetic photometry using 2M2139's NIR. spectrum 
IBurgasser et al.l I2006D f r om the SpcX Prism Library. Note that 
the closely agreeing 2003 and 2004 epochs have been averaged. 
"= R,angcs span the entire set of WIRC observations (i.e. all epochs) 
presented here. 

"* R.anges correspond only to epochs where near-simultaneous JHKs 

photometry is available. 

"= Derived from the empirical relation of IStephens et al.l i2009D 
Derived from the empirical relation of lMarocco et al.l I I2010I) 

band peak, while having a much smaller effect on the H 
and K band SEDs. However, even this composite match 
fails to reproduce the relative strength of CH4 absorption 
in the J, H and Kg bands, leading the authors to con- 
clude that if a binary, the components "may themselves 
have unusual properties." Given the unique variability 
of 2M2139 reported in the present work, its peculiar NIR 
spectrum may be the result of an uncommonly heteroge- 
neous atmosphere, rather than binarity. 

2MJ2139 was observed on 2006 June 23 using 
HST/NICMOS with the NICl camera (43 mas pixel"!) 
and the F170M and FllOW filters (HST program ID 
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Figure 7. J (black) H (red), and Ks (blue) light curves span- 
ning 01 Aug 2009 to 08 Nov 2009 shown on a common flux scale 
where unity corresponds to Ja/xOi H, and Kg magnitudes of 14.75, 
14.11, and 13.59 respectively. The data have been phased to an 
over-constrained period of 7.723 hr (see text for the explanation). 
For clarity, the data have been binned by factors of 3 (WIRC J- 
band sequences), 7 (CPAPIR J- band sequence in bottom panel), 
and 5 (H and Kg sequences). The interleaved JHKs sequences are 
binned as in figure U The best-fitting two-spot model from figure 
[5] is overplotted (solid grey lines) to highlight deviations from the 
original light curve shape as a function of time. The same model, 
with the amplitude scaled by 55%, is also overplotted on the lower 
amplitude sequences as a visual aid. 



10143). We analyzed these data to put constraints on 
the binarity of this brown dwarf. For point-spread func- 
tion (PSF) fitting and subtraction, we used both a model 
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Figure 8. Light Curves for 2M2139 spanning 21 Sep 2009 to 01 
Oct 2009 phased to over-constrained (see text for explanation) pe- 
riods of 7.723 hr (top) and 15.446 hr (bottom). Black, red and blue 
points represent J, H, and Ks band data respectively. Different 
symbols of the same color differentiate different epochs. The data 
have been binned for clarity. Unity corresponds to Jmko> H> ^^'^ 
Ks magnitudes of 14.75, 14.11, and 13.59 respectively. 

PSF generated with the TinyTin:0 software (jKrist 199^ 
as well as the image of another brown dwarf (2MJ0257- 
3105) observed with the same settings only two days 
later as part of the same program. The latter refer- 
ence PSF provided a noticeably better fit to the image 
of 2MJ2139+02 than the former, although in both cases 
the subtraction left no significant residuals. Thus the 
data provide no indication of binarity. By introducing 
fake binary companions of various contrasts and sepa- 
rations in the images and repeating the reference PSF 
fitting, we can rule out the presence of a binary compan- 
ion with contrasts of 0, 1, 2, and 3 mag at separations 
larger than 0.055", 0.065", 0.08"and 0.13", respectively. 
These limits correspond to physical separations between 
0.66-1.56 AU at a distance of 12 pc. 

Given the apparent proximity of this source, future 
parallax measurements should provide more stringent 
constraints on its distance, absolute magnitude, and po- 
tential binarity. 

4.2. Model Fits to the NIR SED 

Detailed atmosphere models including the effects of 
dust condensation and settling can be used to explore 
atmospheric properties of field BDs including effective 
temperature, gravity, metallicity, vertical mixing (via 
departures from chemical equilibr ium), and condensate 
clouds (e.g [Gush ing ct al,, _200_^ [Stephens etlUI 120091 : 
IWitte et al.ll201]l ). However, certain issues with atmo- 
sphere models are well documented. For instance, stud- 
ies of benchmark objects — those with known masses (e.g. 

^ |http: //www. stscl ■ edu/hst/observatory/f ocus/TinyTlm| 
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iDuDuv et all 1201 1[ ). and/or ages ()Leggett et al.ll2008D — 
have demonstrated that effective temperatures deter- 
mined from model atmosphere fitting can differ by up to 
a few hundred degrees compared to those derived from 
evolutionary models. Second, low-gravity L dwarfs (Fa- 
hert y et al. 2011, submitted) and directly ima ged planets 
(e.g. iBowler et all 120101: iSkemer et all 120 lit) have been 
found to be less luminous than "normal" field BDs of a 
given spectral type. This latter observation is surpris- 
ing as lower gravity objects are expected to be larger in 
radius, and hence brighter. The ability of low gravity at- 
mospheres to retain a thicker and higher cloud layer, high 
metallicity, an d/or non-equilibrium chemistry may be re- 
sponsible (e.g. iBarman et ani201lD . and highlights that 
significant interdependencies between effective tempera- 
ture, condensate properties, surface gravity, and other 
secondary parameters may exist. 

In addition to known problems, it is unclear whether 
standard ID models are appropriate for 2M2139 as 
it's large variability indicates a heterogeneous surface. 
Nonetheless, we have performed fits of ID model atmo- 
spheres to 2M2139 's NIR spectrum in order to determine 
the general properties required to reproduce its spectrum 
in a manner consistent with other studies. Also, in the 
following section we model 2M2139 's variability and NIR 
spectrum simultaneously using linear combinations of ID 
cloudy and clear model atmospheres. Thus, the single 
spectral fits performed here will serve as a reference to 
which hybrid cloudy/clear models can be compared. 

4.2.1. The Atmosphere Models 

The model spectra employed here have been calculated 
for solar metallicity and chemical equilibrium, for sim- 
plicity allowing only effective temperature (T^ff), sur- 
face gravity (logg) and cloud sedimentation efficiency 
(/sed) to vary. The fsed parameter describes the efficiency 
of condensate sedimentation wit hin the atmosphere, 
accord ing to the cloud model of lAckerman fc MarlevI 
(|2001| ). which has been successfully used to repro- 
duce observations of a wide range of clou dy BDs (e.g. 
iCushing et all 120081: [Stephens et aD l2009| ) as well as 
Jupiter. Large values of /sod correspond to more effi- 
cient particle growth, resulting in larger particle sizes 
that more efficiently rain out of the atmosphere and con- 
sequently lower opacity cloud decks. Our model grid was 
calculated for 500 K< T^jf <2000 K in steps of 100 K, 
log5 ={4.5,5.0}, and /sod={l,2,3,4,nc}, where 'nc' de- 
notes a cloudless model. Although the models were cal- 
culated with 100 K spacing, we have linearly interpolated 
between them to obtain a finer sampling in T^ff of 25 K. 

4.2.2. Fitting procedure 

For each model in our grid we computed a weighted 
statistic, foll owing a modified ver sion of the procedure 
described by ICushing et al.l (|2008() , and given by 

G = ^u;(AO f Pj^^) -^M{X,;T,ff, log gJsed) Y 

where F{Xi) is the observed spectrum, and 
M{Xi;Teff,\ogg, fsed) is the model, scaled by a 
constant factor c = {R/d)^ that minimizes G, where R 



is the BD radius and d is its distance. All model spectra 
were smoothed to replicate the average resolution 
of 2M2139's spectrum over the wavelength region of 
interest (R^lOO), via convolution with a Gaussian. Due 
to unequal spacing of the data, the fit is weighted by 
w{Xi) — AAi, where AA^ is the width of wavelength bin 
i. Fitting was performed over the wavelength intervals 
of 1.0-1.35 /J,m, 1.42-1.8 ^m, and 1.95-2.35 /im, avoiding 
regions of high telluric water absorption. Due to 
uncertainties in molecular band strengths and opacities, 
as well as the complexity of condensate cloud physics, 
the models are generally incomplete and can disagree 
with the observations at a level significantly greater 
than the measurement uncertainties associated with the 
data. Thus we have chosen to set all = 1, so as not 
to bias a particular wavelength regime based on their 
measurement uncertainties. This is equivalent to includ- 
ing a constant uncertainty term for each wavelength bin 
in order to account for model incompleteness. Since 
this constant term is much larger than the measurement 
uncertainty for any given wavelength bin included in the 
fit, the Ui become essentially constant. 

In figure |9| we plot the G statistic (normalized by its 
absolute minimum value across all models tested, Gmm) 
for various values of /scd and log 5 = {4.5,5.0} as a 
function of effective temperature. We find the best fit 
corresponds to logg = 4.5, /scd=3 and Te//=1400 K. 
While an effective temperature of 1400 K is consis- 
tent with the range of L/T tra nsitio n temperatures 
found by iGolimowski et al.l (|2004[ ) and [Stephens et al.l 
(|2009[ ). a surface gravity of log (7=4. 5 would then imply 
a rather young age of 100 Myr and mass of ~20 Mj^p 
(jSaumon fc MarlevI 2008). It is interesting to note that 
the K band portion of 2M2139's SED has been reported 
as a best match to that of the low gravi ty directly imaged 
planet IIR8799 b (jBarman et al.ll2011| l. However, model 
fits to the low resolution NIR spectra of field BDs do not 
strongly constrain surface gravity, and can yield values 
of log (7 differing by up to 1.0 dcx when different wave- 
length regimes are used for fitting (Gushing et al. 2008). 
In addition, according to IGushing et al.l (2008) a value 
of log (7 = 4.5 derived from spectral fitting is not par- 
ticularly unusual for a field BD, even when evolutionary 
sequences predict higher surface gravities (log (7 > 5) for 
the same objects. Thus, it would be premature to con- 
clude anything about 2M2139 's age and gravity based 
on model fits alone, especially since there is no addi- 
tional evidence for youth. With this in mind, it may be 
more likely that interdependencies between surface grav- 
ity and other physical parameters such as cloud thickness 
and/or metallicity (despite being assigned independent 
model parameters) lead to the latter being mimicked by 
the former in spectral models. For the purpose of this 
paper we continue to use the log (7 = 4.5 models to de- 
scribe 2M2139 's atmosphere, but caution that the phys- 
ical interpretation of the gravity parameter is unclear. 
The analysis presented here and in subsequent sections 
is not strongly affected by this choice of model surface 
gravity, other than providing an optimal fit to 2M2139 's 
NIR SED. Specifically, trends pertaining to photomet- 
ric variability in the following sections are qualitatively 
similar for different surface gravities. Furthermore, from 
figure |9l Teff appears roughly independent of log g, but 
depends strongly on fsed- Thus, while log (7 is poorly 
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Figure 9. Left panel: The weighted statistic G (normalized to its overall minimum value, Gmin) for a series of different /gedi 'ogg, 
and T^ff. Right panel: The model spectrum corresponding to the minimum G statistic (green line) is shown ovorplotted on the data (black 
line). Regions of high telluric absorption which have been omitted from the fit arc shown in grey. 



constrained by the models, it is only of secondary impor- 
tance to our analysis in comparison to the cloud thickness 
parameter, /sod- 

5. MODELING VARIABILITY DUE TO HETEROGENEOUS 
SURFACE FEATURES 

While wc discuss alternative explanations for 2M2139's 
variability in the discussion section, here we construct a 
simple model for the observed variability, assuming that 
heterogeneous surface features are responsible. We envi- 
sion a scenario where "spots" or cloud features remain 
approximately static in the BD's rotating frame over a 
single rotation, and variability arises due to rotational 
modulation. We consider a simple model where the BD 
surface is composed of two types of regions differing in 
temperature and/or cloud properties, possessing surface 
fluxes Ti and J^2- At a given snapshot in time the total 
flux from the BD, is given by a linear combination 
of J^i and J-2 weighted by their relative filling fractions 
over the BD's visible disc. The peak-to-peak amplitude 
of variability that an observer would detect in a given 
bandpass due to a change in these filling factors can be 
expressed as 



(1 



Aa) J"i + (a + Aa) J"2 - (1 - a) J"i - aj^2 



0.5[(1- 



■ Aa)J-i + {a + Aa)T2 + (1 ■ 
Aa 



(2) 



(3) 



where A = /F is the change in flux divided by the 
mid-brightness flux. The parameter a is the minimum 
filling factor of the F2 regions, Aa is the change in filling 
factor, and AJ^ = ^2 ^ ^i- The parameter a=a-|-0.5Aa 
corresponds to the filling factor of the ^"2 regions at mid- 
brightness. 

Synthetic photon fluxes T2 and J-i for use in equa- 
tion [3] were computed from the ID model spectra in- 
tegrated over the WIRC system plus filter transmis- 
sion curves. For t his pu rpose we used the models of 
iSaumon fc MarlevI ()2008[ ) described above with solar 
metallicity, log 5 = 4.5, and a range of values in fseA 



and Te/ / . We also experimented with different log g val- 
ues, but found them to be a poor match to 2M2139 's 
NIR spectrum, and hence only present the optimal case 
where log g=4.5 here. In general, J^i = J^i[T'i, /scdi] and 
F2 = J-2 [2^1 + AT, /sod2] so that each region is character- 
ized by a distinct effective temperature and cloud sedi- 
mentation efficiency. As a matter of convention we will 
always choose /sodi < /sod2, such that the Fi regions 
have higher condensate opacity. The parameter Ti rep- 
resents the effective temperature of the Fx regions, while 
Ti -I- AT is the effective temperature of the less cloudy 
F2 regions. For a given combination of Ti, AT, /sodi, 
/sod2, a and Aa, fluxes can then be determined from 
ID model atmosphere grids, and the resultant variabil- 
ity computed via equation[3l We caution that, in general, 
ID models with different effective temperatures and ver- 
tical distributions of dust condensates have different un- 
derlying pressure-temperature profllcs, and thus interpo- 
lations between them cannot be made in a self-consistent 
way. Nonetheless, in the absence of 3D models, it is in- 
structive to use ID models as a guide, understanding 
their drawbacks. Limitations of our modeling approach 
are discussed further in section 

Using equation [3] we computed AkJAj and Ah/Aj 
as a function of AT for a variety of combinations of /sodi 
and /sed2- The various combinations can be broadly di- 
vided into two main cases, which we specify for later 
reference: 

Case A — Case A corresponds to an atmosphere with 
heterogeneous cloud features, where /sodi 7^ /sed2- This 
case could describe an atmosphere composed of clouds 
and clearings, or more generally regions of differing 
condensate opacities. This picture of a brown dwarf 
with spatially variable cloudiness is motivated by the 
cloud fragmentation hypothesis of the L/T tran sition 
(|Ackerman fc MarlevI [200TI : iBurgasser et al.ll2002| ). dis- 
cussed in the introduction. 

Case B — Case B corresponds to a uniformly cloudy 
atmosphere where /sedi = /sed2, with heterogeneities in 
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temperature only. This case may approximate the pres- 
ence of cool or hot magnetically induced spots. We have 
not modeled hot or cool spots in a uniformly clear at- 
mosphere because the red J — Kg color of 2M2139 is 
incompatible with cloud-frce models and precludes this 
scenario. For this case we have chosen to consider only 
/sodi = ./scd2 = 3 based on the best-fitting ID model 
atmosphere for 2M2139 (figure [H). 

For each pairing of model grids specified by 
/sedi={l,2,3} and /sod2 = {3,4, nc} we used equation[3] 
to compute {\)AkJAj and A^/Aj as a function of AT, 
and (ii)thc corresponding change in filling factor, Aa re- 
quired to produce the maximum observed variability of 
Aj=0.26. Results for a representative selection of these 
pairings arc shown in figure [10] for values of ri=1100 K 
and 1400 K; a reasonable range for 2M2139. Results are 
plotted for both a = 0.2 and a = 0.6 in order to demon- 
strate a weak dependence of our result on the minimum 
filling factor. Note that the overall effective temperature 
can be much greater than Ti when AT > 0, as warm 
regions quickly dominate the total flux. Although it is 
not clear a priori what appropriate values of Ti should 
be, 1400 K represents a reasonable upper limit when 
AT > 0, since it coincides with both our model fitting 
in section 14.21 and the upper limit of the L/ T transition 
temperature found by [Stephens et ahl (|2009[ ). The lower 
value of Ti=1100 K is shown to be appropriate in the 
following section. 

For comparison, figure [TT] compares results for case A 
from different model atm osphere groups , including the 
cloudy and clear models of lBurrows et al.l (I2006D and th e 
SETTL and COND models of lAhMd et al .1 (200lL[200l . 

Our figures IIOIII II reco ver trends previously discussed 
by lArtigau et al. ~ |2009l ) . Namely, cool or hot spots 



within a cloudy atmosphere (case B) produce ampli- 
tude ratios Ak,/Aj exclusively > 1, whereas warm 
regions of low-condensate opacity (case A) yield < 
AkJAj < 1 for AT > 100 K . The observed ampli- 
tude ratios for 2M2139 are positive and consistently <1 
and therefore inconsistent with the presence of magnetic 
spots. Rather, our modeling supports an interpretation 
of patchy clouds, wherein regions of lower condensate 
opacity are warmer. 

As seen in figures [TUlfTTl a variety of model combi- 
nations for case A are able to reproduce observations, 
with differing requirements for AT and Aa. In general, 
models where the difference between /sedi and /sod2 is 
small require lower temperature contrasts but sightly 
larger changes in filling factor to model the observa- 
tions. Increasing the parameter a yields large ampli- 
tude ratios for a given temperature contrast AT, but 
also requires a larger Aa. In all cases temperature con- 
trasts > 150 K are required to model the observed 
amplitude ratios, with some model combinations re- 
quiring temperature contrasts in excess of 400 K. The 
SETT L/C OND and cloud y /clear models of lAllard et all 
()2003[) and iBurrows et al.l ()2006[ ) are equally able to re- 
produce the photometric variations, requiring somewhat 
lower temperature contrasts when Ti < 1200 K. While 
the variety of models and parameters capable of repro- 
ducing the photometric data appears large, in the follow- 
ing section we use 2M2139's NIR spectrum as an addi- 
tional constraint to significantly reduce the span of good 



solutions. 

5.1. A specific variability model for 2M2139 

The modeling of photometric variations alone, as done 
in the previous section, ignores additional information 
about 2M2139's spectral energy distribution that can fur- 
ther constrain the nature of the v ariability. In addition , 
it has been previously noted fe.g. IBurrows et al]|2006() . 
that the same models capable of reproducing observed 
changes in broadband colors do not necessarily repro- 
duce the NIR spectra of L/T transition dwarfs. 

For combinations of /sodi = {1,2,3} and /sod2 = 
{3,4, nc} and logp = 4.5, we searched for combinations 
of Ti , AT, a, and Aa capable of reproducing both the NIR 
SED and observed variations of 2M2139. We have re- 
stricted ourselves to the region of parameter space where 
AT > corresponding to the upper right quadrant in 
figure [TOl We systematically examined a grid of hybrid 
models with Ti = 800 - 1600 K and AT = - 700 K in 
steps of 25 K. Our warmest model has Tf.ff = 2100 K, 
and therefore our grid is incomplete for regions where 
Ti + AT > 2100 K. We note that aU original spectra were 
computed with 100 K spacing, and we have linearly in- 
terpolated between models to obtain our finer sampling. 
At each grid point: 

1. A parameter a (corresponding to the mid- 
brightness filling factor, a-|-0.5Aa) was uniquely de- 
termined by fitting a hybrid model given by 



M(A) = (l-a)J-i[A;Ti,/sedi] 

-FaJ-2[A;Ti + AT,/,ed2] 



(4) 



to 2M2139 's NIR spectrum, using a procedure 
identical to the one used for spectral fitting in sec- 
tion 14.21 (sec equation [1]). We note that in the 
above equation we have used J-i and J-2 to rep- 
resent model fluxes as a function of wavelength, in- 
stead of integrated over a given bandpass as before. 
The minimum weighted statistic, G, associated 
with each grid point was recorded. 

2. We solved for a and Aa by inputing a = a — 0.5Aa 
and Aj=0.26 into equation [3l where we have cho- 
sen to set Aj to the maximum observed variabil- 
ityQ In some cases Aa derived from the best spec- 
tral flt was unphysical (e.g. a -I- Aa > 1), and in 
these cases we systematically varied our fit param- 
eter a until we found the lowest-G fit that also 
satisfied the physical condition that a + Aa < 1. 

3. We computed model amplitude ratios AkJAj and 
Ah/Aj by inputting parameters a, Ti[Ti], and 
7^2 [Ti + AT] into equation [3l A statistic for 
the amplitude ratios given by Xamp = [{^h/Aj — 
x)^/i7^ + i^Ks/^J — y)^/'^l] was recorded, where 
X and y are the observed amplitude ratios, and ax 
and cFy their uncertainties (sec tabic [3|). 

When a is fixed, Aa is linear in Aj, and therefore the Aa 
corresponding to smaller overall amplitudes can be easily inferred. 
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Figure 10. Top: Model amplitude ratios AkJAj for a BD with a heterogeneous surface as a function of the temperature 
contrast between surface elements, AT. Amplitude ratios were computed from equation [3l using a variety of cloudy and clear 
model fluxes ( J^i [/scdi , 21] and T2 [/sed2 , T2 + AT] ) from Saunion & Marlcv (200^) , indicated by color in the plot legend. Here we 
have used the shorthand of fni/fn2 to denote /scdi ~ "-i//sGd2 ~ n2, and 'nc' for the cloud- free model. For all cases curves are 
plotted for constant ri = 1100 K (solid lines) and Ti = 1400 K (dotted lines). Colored lines correspond to case A (heterogeneous 
clouds) while the black lines correspond to case B (cool/hot spots). The left and right panels demonstrate a weak dependence 
on the initial filling fraction of the J-2 regions, with a— 0.2 (left) and a=0.6 (right). Middle: For each curve in the above panel, 
the corresponding change in filling factor, Aa, required to produce the maximum observed variability for 2M2139 of Aj =0.26, 
is shown. Bottom: Same as the top panels but for Ah/Aj {Aa2e remains unchanged from the upper panel). 
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Figure 11. Same as figure 1101 with a = 0.2, and including 
comparisons t o the clou dy /clear and SETTL/CO ND models of 
IBurrows et"all pOOBl ) and lAllard eFall l(200ll . \2003t ) ( sec the main 
text for more detail). For each set of models the three curves cor- 
respond to fixed values of Ti=1100 K, 1200 K, and 1400 K. 

According to the above procedure (repeated separately 
for various combinations of /scdi and /scd2), G and ximp 
were determined for every grid point in the Ti-AT plane, 
allowing us to identify regions of parameter space capable 
of reproducing both the NIR spectrum and broadband 
variability of 2M2139. 

We obtained best-fit parameters by determining where 
Xamp+4 G/Gsingic was a minimum on our grid. Here, 
Gsingio corresponds to the best-fitting single (non-hybrid) 
spectral model found in section |4?2] We note that the G 
values have iV— 3 degrees of freedom (dof; where N = 304 
is the number of data points used for fitting) and have 
been scaled such that Gsingio = dof. If we were to mini- 
mize Xamp+^i the best-fit parameters would be strongly 
dominated by the spectral fit, which uses significantly 
more data points (however, due to model incomplete- 
ness and highly correlated residuals this increase in data 
points does not result in proportionately more precise 
parameters). Instead, we have used our prior knowledge 
of the observed spectrum to impose a more broad con- 
straint, allowing a wide range of "credible" spectral tem- 
plates and eliminating obvious mismatches. In practice, 
we found that scaling G by 4 / Ggingie assigns a reasonable 
weight to the spectral fit (e.g., see figures [T^ and [T^ . 

We found only a few combinations of /scdi and /sGd2 
that were capable of reproducing the photometric vari- 
ations, while at the same time providing a reasonable 



spectral match. Those where cool, high condensate 
opacity regions are represented by /sedi = {1,2} and 
warm, lower condensate opacity regions are represented 
by /sed2 = {3, 4} generally provide good fits to the data, 
with the /sodi = l//sod2 = 3 combination providing the 
best fit. These four model scenarios, and corresponding 
best-fit parameters, are shown in figure [T^ fitted to the 
simultaneous JHKg photometric data from 30 Sep 2009. 
Shown are contours for AG/Gsingie and Xamp tti^ ^i" 
AT plane. The hybrid model spectrum corresponding 
to the best-fit parameters, as well as the resultant vari- 
ability as a function of wavelength are also shown. We 
have chosen to show fits to the Ah /Aj and AkJAj am- 
plitude ratios from 30 Sep 2009 as the amplitude ratios 
measured on this date are approximately intermediate 
to those measured in the other two epochs for which we 
have simultaneous JHKg light curves. We have opted 
not to average the amplitude ratios over all epochs since 
the measurements are incompatible at the 2a level and 
may represent different cloud configurations. In partic- 
ular, the near-unity amplitude ratios from 26 Sep 2009 
are challenging to fit, and strongly bias the averaged val- 
ues. Rather, we provide best-fit contours and spectra 
corresponding to the other two epochs in Appendix |B] 

Depending on the combination of /scdi and /sod2 we 
infer temperature contrasts of AT=175-425K between 
thick and thin cloud regions, and changes in filling factor 
of Aa=0.12-0.39 (see figurelHl). 

Rather than an atmosphere composed of clouds and 
clearings, our observations are best reproduced by one 
possessing regions of differing cloud opacity. For the 
best-fitting pairing of /sodi = l//sed2 = 3 shown in figure 
1121 (and similarly for the /scdi = 2//sed2 = 3 combina- 
tion) the data can be described by an atmosphere com- 
posed primarily of /scd2 — 3 regions (67% -93% coverage 
over the visible disc as a function of phase) with cooler, 
thicker cloud patches (/scdi = 1) covering the remainder 
of the surface (7%-33%). The thick cloud patches have 
a temperature of ri = 1125 K, while the thinner cloud 
patches are AT = 250 K warmer. Physically, this model 
could indicate condensate clouds that form in at least 
two distinct layers at different depths in the atmosphere. 
Such layers could arise d ue to compositional stra tifi- 
cation of condensates (e.g. iLodders fc Feglevl I2006D . or 
from complex atmospheric dynamics. Long-lived cy- 
clonic storm systems and transient high altitude clouds 
are seen in the atmospheres of the Solar System gas gi- 
ants fe.g. lSmith et al.lll989tlVasavada et al.iil998D . and it 
is possible that similar features exist in the atmosphere 
of 2M2139. Since small spatial scales will tend to be 
averaged out in the disc-integrated light, the large am- 
plitude of our observations (and the corresponding large 
azimuthal asymmetry in modeled cloud coverage) is most 
consistent with the presence of sizable cloud features. If 
the observed variations were attributed to a single high- 
altitude storm, it would have to occupy 26% of the visible 
disc in order to produce Aj =0.26. 

Alternatively, for the next-best-fitting /scdi = 
l//sod2 = 4 model also shown in figure [T^ (and simi- 
larly for the /scdi — 2//sod2 — 4 combination) the data 
can be described by an atmosphere composed primarily 
of cloudy /scdi =1 regions with Ti=1025K (67%-55% 
coverage as a function of time) interspersed with warmer 
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Figure 12. Simultaneous fits to both the NIR spectrum (black line) and photometric variability (black data points) of 2M2139 for the 
four model combinations where /aedl = {Ij 2} and /sod2 = {3, 4} (see text for explanation). Right: Contours in the Ti-AT plane for Axamp 
(blue) and for 4: AG/Gsingie (red), corresponding to levels of {2.3,6.2, 11.8}, where the A indicates differences in G and Xamp from their 
minimum values on the grid. A red cross marks the combination of Ti and AT (for which a and Aa are also uniquely determined) where 
Xamp+4 G/Gsingie is a minimum. The shaded map shows the exponential of — (Xamp +4 G/Gsingie)/2 with square-root scaling to emphasize 
the tails of the distribution. Left: The hybrid model spectrum at mid-light (red solid line), corresponding to the best- fit parameters marked 
by the red cross on the right. The model is a linear combination of cloudy (cyan line) and less-cloudy (green line) surfaces, where the 
filling fraction of less-cloudy regions varies by an amount Aa. Dotted red lines show the hybrid model at maximum and minimum light, 
assuming Aj=0.26, and normalized over the 1.23 to 1.25 fim range. The bottom left panel shows the relative variations that would result 
for the best-fiting model as a function of wavelength, normalized to the J band. The observed amplitude ratios and their errors are over 
plotted for comparison, and Xamp for t^^e best fit is specified. 
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regions of lower condensate opacity (/scd2 = 4). In this 
case a much higher temperature contrast of AT = 425 K , 
and a correspondingly smaller change in filling factor af- 
fecting only 12% of the visible disc are required. This 
scenario is closer to the picture of an atmosphere cov- 
ered by cloudy and clear regions that might naively be 
expected based on theories of cloud fragmentation at the 
L/T transition. However, it is notable that we are un- 
able to find solutions where /sed2='nc' to represent the 
data, indicating that the modeled "clearings" must retain 
some degree of condensate opacity (i.e. /scd2 = 4). For 
an atmosphere with a same T-P profile throughout and 
where temperature decreases monotonically as a function 
of altitude (as is expected for isolated BDs), this large 
temperature contrast is an indication that thick cloud 
patches extended to higher, cooler regions of the atmo- 
sphere. 

For the best-fitting hybrid spectra, the normalizing 
constant c = [R/d) m Eq. [Tl along with a BD ra- 
dius of i?=0.13 obtained f rom the evolution mod- 
els of iSaumon fc MarlevI ()2008l) yield spectroscopic dis- 
tances ranging from 15.7-18.5 pc, which are larger than 
the 11.6 pc distance deri ved earlier fr o m th e empirical 
K-hand relationship of Ma rocco et al.l (|2010( ). If we in- 
stead assume a higher gravity of log g = 5.0 at odds with 
the best-fitting spectral models, we find R — 0.1 Rq and 
obtain distances of 13.8-14.1 pc, which are closer to the 
empirical estimate. 

We also attempted the above analysis using the 
SETTL/COND an d clou dy/clear model combinations of 
lAllard et all (|2003[ ) and iBurrows et al.l (|2006[) . however 
we were unsuccessful in simultaneously reproducing the 
NIR SED and photometric variations of 2M2139 using 
these models. This is not surprising given that we find 
no good solutions with /scd2 = 'nc' (truly condensate-fr ee 
regions) using the models of ISaumon &: MarlevI ()2008l ). 

Finally, we should note that in finding the best hy- 
brid models, we have assumed that our target's observed 
spectrum corresponds to its mid-brightness state. This is 
a reasonable approximation, as synthetic colors derived 
from 2M2139's NIR spectrum are consistent with J — Kg 
and J — H colors measured for this target at the time 
of the WIRC observations. In figure [T^ we display the 
model maximum and minimum light spectra in addition 
to the best-fitting mid-brightness spectrum, all normal- 
ized to the J band, in order to demonstrate that differ- 
ences between the three are minimal. Even when the 
absolute variability amplitudes are large and differ as a 
function of wavelength (Aj=0.26, AkJAj = 0.45-0.83), 
changes to the shape of a normalized spectrum are barely 
perceptible by eye. Accordingly, the assumption that the 
observed spectrum corresponds to mid-brightness does 
not significantly impact our results. In contrast, the ob- 
vious spectral mismatch that results from fitting models 
to photometric variations alone, without using a spectral 
constraint, is illustrated in figure [T^ for SIMP 0136. 

5.2. Comparison to SIMP0136 

Here we apply our variability model to SIMP0136, 
the oth er known example of a quasi-periodic variable T- 
dwarf. lArtigau et al.l ()2009D modeled the variations of 
SIMP0136 {AkJAj ~ 0.48 and Aj - 0.05 mag) us- 
ing ID cloudy and clear model atmospheres and an in- 
terpolation scheme equivalent to our equation |3l While 



we ha. v e used primarily the models of ISaumon fc MarlevI 
(|2008f) . lArtigau et aLl pOOgj ) used combinations of cloudy 
and clear models from 3 different group s: the SETTL 
and COND models of lAllard et all (l200ll IpOOl and the 
cloudy and cl ear mo dels of Burrows ct al.' (20Q3), as well 
as those from lTsuii fc Nak aiima (2003) . They found that 
a temperature contrast of AT ~ 100 K between clear and 
cloudy regions and Aa ^0.15 can reproduce the J and Kg 
variability of SIMP0136, with good agreement between 
all model sets. However, we find that the model spectra 
corresponding to this result (Ti=1100 K, AT=100 K, 
a=0.25) provide a poor match to SIMP0136's NIR spec- 
trum (figure [13]) , demonstrating that models capable of 
reproducing changes in broadband colors do not neces- 
sarily yield good spectral matches. 

In order to make a fair comparison between 
2M2139 and SIMP0136 wc have repeated our above anal- 
ysis for SIMP01 36, relying on t h e amp litudes of variabil- 



ity reported by lArtigau et all (12009 ) and SIMP0136's 
NIR spectrum from iBurgasser et al.l (|2006f ). We find 



that the variability and NIR spectrum of SIMP0136 can 
be well-reproduced by a wide range models, including 
those where /sedi = {1,2,3} and /sed2 =4, as well as 
those where /sedi = {3,4} and /sed2 = nc. These solu- 
tions require temperature contrasts ranging from ~75 K 
for the /scdi = 3//scd2 =4 pairing to ^--^275 K for the 
/scdi = 3//sod2 =nc pairing, and require changes in fill- 
ing factor affecting 3% to 15% of the visible disc. We also 
repeat ed our analy s is usi ng the SETTL/COND mod- 
els of lAllard et all (I2003D and cloudy/clear models of 



IBurrows et al.l ( 20061 ). For these latter models we found 
poor overlap between parameters capable of simultane- 
ously reproducing the observed photometric variations 
and SIMP0136's NIR spectrum, with optimal tempera- 
ture contrasts occurring for AT=500 K and 300 K re- 
spectively, and significantly poorer spectral fits. A se- 
lection of best-fitting hybrid spectra for SIMP0136 from 
different model groups is shown in figure [131 

In general, we find that ou r result s agree qualitatively 
with those of lArtigau et al.l (j2009() . However, for the 
sake of comparison, we find that temperature contrasts 
required to model the observed AkJAj ~ 0.48 for 
SIMP 0136 are at least 2-3 times higher (and changes in 
filling factor smaller) than previously inferred if a spec- 
tral constraint is imposed. 

5.3. The need for self-consistent 3D model atmospheres 

Due to the lack of self-consistent 3D models our above 
modeling approach is necessarily data-driven as opposed 
to physically motivated. The drawback is that combi- 
nations of independent ID atmosphere models may pos- 
sess incompatible T-P profiles, and we are therefore not 
guaranteed to arrive at physi cally plausible con figura- 
tions (see the discussion within iMar lev et al.ll2010( ). This 
leads us to question whether the high temperature con- 
trasts derived above, and in particular AT > 300 K in 
several cases, are physically reasonable. 

Figure[T4lshows the temperature of the t = 2/3 surface 
as a function of wavelength for a model with /sod = 2, 
log (7 = 4.5 and Te//=1000 K, with and without conden- 
sate opacity. Both models share the same T-P profile, 
but in the latter condensate opacity has been switched 
off. The cloud- free r = 2/3 surface occurs at higher tem- 
peratures, indicating that photons emerge from deeper 



18 



Radigan et al. 



1 aL SIMF^0136 

■ Allard etal. SETTL(1100 K)/COND(1200 Kl, a=0.25 

- Burrows et al. cloudy(1 100 K)/clear(1 200 K), a=0.25 

1.2 - 




X ()im) 



1 4L SIMF^0136 

- SETTL(1350 K)/COND(1850K),a=0.13,Aa=0.02 

- Burrows et al. cloudy(1275 K)/clear(1575 K), o[=0.22, A a=0.03 

, n Marley etal. f3(1250 K)/nc(1525 K), a=0.12, A a=0.03 

^■'^ Marleyetal.f1(1100)/f4(1350 K),ci=0.89,Aa=0.05 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 

X (|j,m) 



Figure 13. The NIR spectrum of SIMP0136 IIBurgasser et ani2006l ). with hybrid model spectra ( line ar combinations of cloudy and clear 
models) used to model its photometric variability overplotted (similar to the left panel of figure [12}. The filling fractions and effective 
temperatures of the models used to form each hybrid spectrum are shown in the legend. Left: Hybrid model spectra obtained by setting 
Ti + AT=1200 K a priori, and modeling only photometric variations as done bv lArtigau et aLl l|200g| ') (see text for more detail). Right: 
Hybrid model spectra corresponding to our simultaneous modeling of the photometric variations and NIR spectrum of SIMP0136. Red 
and green lines correspond to the same model sets as in the left panel, but with filling fractions and temperature contrasts that better 
re produce SIMP0136 's NIR spectrum as well as its photometric variability of Aj^^/Aj = 0.48. Blue and pink lines correspond to models 
of jSaumon fc MarlevI II2008I ') where f^adi = l//sod2 = 4 and f^^di = /4//scd2 ='nc' respectively. 



atmospheric layers in the absence of condensate opac- 
ity. The temperature contrast between the cloudy and 
cloud-free "photospheres" is as high as 600 K in the J 
band, and therefore the large temperature contrasts re- 
quired to model our observations are not unreasonable. 
However, the relative contrasts in the J and Ks bands 
appear problematic; in the Kg band the contrast be- 
tween the cloudy and cloud-free r = 2/3 surface has 
dropped to only ~100 K. Thus, in the Ks band, one 
does not look much deeper into the atmosphere even if 
cloud opacity is removed. This is mainly due to an in- 
creased gaseous opacity from collision induced absorption 
by H2 molecules (CIA) at these wavelengths, causing the 
gaseous (cloud-free) photosphere to occur at higher alti- 
tudes, and hence cooler temperatures. This structure of 
the cloudy model poses a problem if we wish to reproduce 
the large observed amplitude ratios of Ak^/Aj > 0.45, 
as it requires the Kg band flux within clear regions to 
be larger than allowed by the T-P profile of the cloudy 
model. Of course, this assumes that the T-P structure 
of the upper atmosphere is identical within cloudy and 
clear regions, which need not be the case. 

The above results suggest that conclusions drawn from 
linear combinations of ID atmosphere models must be 
viewed with caution, and highlights the need for self- 
consistent 3D modeling, including the effects of atmo- 
spheric circula tion. Recent progre ss in this direction has 
been made bv lMarlev et al.l (|2010f ). who have calculated 
a set of self-consistent patchy cloud models of the L/T 
transition, where cloudy and clear regions share a com- 
mon T-P profile. These model s conti nue to use the cloud 
model of lAckerman fc MarlevI (j2001l ) to describe cloudy 
regions and contain an additional parameter h, which 
specifies the fractional coverage of "holes" , or cloud-free 
regions in the atmosphere. However, these models pre- 
dict negative Ak^/Aj (i.e. anti-correlated amplitudes), 
inconsistent with th ose observed for 2M2139. Since the 
iMarlev et aP ()2010[ ) models assume a common T-P pro- 
file within cloudy and clear regions, this may indicate 
that the T-P profiles diverge in a real atmosphere. Al- 



ternatively, our modeling in the previous section demon- 
strates that the best fits to the data derive from models 
in which the lowest condensate opacity surface elements 
are not truly condensate- free (i.e. /sed = 4), which may 
also account for p art of the disc r epanc y between the data 
and models from IMarlev et all ()20 10f). 

There arc several possible explanations for the mis- 
match between self-consistent models and the Kg band 
data. A larger than expected contrast between cloudy 
and clear regions in the Kg band could suggest thick, 
high-altitude dust clouds, or weaker than expected H2 
CIA, possibly the result of high metallicity. Alterna- 
tively, the T-P structure of a partially cloudy atmosphere 
may itself be horizontally inhomogeneous within the pho- 
tosphere, requiring a fully 3D hydrodynamical treatment 
to model. In addition, there is always the possibility 
that different types of cloud heterogeneities than envi- 
sioned here, or a completely different mechanism may 
be responsible for the variability of 2M2139 . Assuming 
that the variability is indeed the result of heterogeneous 
clouds, these observation can be used to guide future 
modeling efforts. 

6. DISCUSSION 

6.1. Origin of the Observed Variability 

While heterogeneous clouds is our preferred explana- 
tion for the observed variability, and the scenario we have 
modeled in depth, it is worth exploring whether other 
mechanisms could be responsible. 

6.2. Clouds versus Magnetic Spots 

The asymmetric and evolving light curve shape of 
2M2139, with periodic behavior on a rotation timescale, 
is the typical signature of surface spots. As noted ear- 
lier, the derived 7.721±0.005 hr periodicity would be 
consist ent with rotation rate s of ultracool dwarfs (^2- 
12 hrs iReiners fc Basril 120081 ) inferred from vsini data, 
whereas a double-peaked light curve with a 15.44 hr pe- 
riod would correspond to a slower rotator (this is not 
inconsistent with v sin i data, as observations are biased 
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Figure 14. Temperature of the t = 2/3 surface as a func- 
tion of wavelength for a model with /s^d = 2, logg = 4.5, and 
Teyy=1000 K, with (black line) and without (red line) condensate 
opacity, using models of lSaumon &: MarlevI lj2008l ). 

against such slow rotators). The nature of the surface 
features responsible is not immediately clear. For stars, 
cool magnetic spots are a common source of variabil- 
ity, related to the suppression of convective heat flux 
within magnetic flux tubes that penetrate the stellar 
photosphere. However, estimates of magnetic Reynolds 
numbers in cool dwa rf atmospheres (jGelino 
iMohantv et al.|[2002l ). predict weak coupling between the 
gas and magnetic field, with magnetic Reynolds num- 
bers of Rm ~ 0.01-10^^ throughout the photospheres of 
late-L and T type BDs. It is unknown to what extent 
cooling in deeper, ionized layers might carry through to 
the photosphere. From an empirical standpoint, only a 
small fraction of late-M dwarfs display variability due to 
magnetic s pots, and typical photometric amplitudes are 
< 5% (e.g.l Terndr up et aT][l99llScholz fc Eisloffe]|[200l 
iScholz etal . 2009). Thus, if magnetic in origin, the NIR 
variability of 2M2139 would be an order of magnitude 
larger than what is typical of late M-dwarfs. 

Alternatively, condensate clouds in rapidly rotating 
BD atmospheres may form discrete cloud features, sim- 
ilar to the banding and cyclonic storms observed on 
Jupiter. At wavelengths where the gaseous opacity is 
low (e.g. in the NIR 1-2 /im regime between water ab- 
sorption bands), the flux emitted from 'holes' or clearings 
in the cloud layer may emerge from atmospheric depths 
hundreds of degrees hotter than the opaque cloud layer 
(e.g. figure [M]). Such large temperature contrasts would 
provide a natural explanation for large-amplitude vari- 
ability in the J and H bands, and somewhat smaller 
variability in Kg- On Jupiter, regions of low condensate 
opacity are visible as bright hot spots at 5 //m, where an 
analogous gaseous opacity window exists. By conduct- 
ing synthetic photomet ry on resolved images of Jupiter, 
iGelino &: MarlevI ()2000[ ) determined that the contrast be- 
tween cloud features and hot spots in Jupiter's atmo- 
sphere may lead to variability as high as 20% at 5/im 
over a single rotation, with variations strongly correlated 
with the position of the Great Red Spot, a high altitude 
anti-cyclone. Thus Jupiter may prove to be the clos- 
est astrophysical analog for 2M2139's unique variability. 



with both warm clearings and a high-altitude storm fea- 
ture contributing to variations. 

Finally, the modeling presented here (SJU figures [TOl 
llip strongly supports the heterogeneous clouds interpre- 
tation — rather than magnetically induced cool or hot 
spots — as it reproduced the observations best. Since 
magnetic structures and clouds would be effectively de- 
coupled, the magnetic spot model would correspond to 
the case B we have investigated (uniformly cloudy at- 
mosphere with spatial variations in T^f / only) , which we 
found cannot explain the observations. This is of course 
consistent with the fact that clouds are known to occur in 
BD atmospheres while magnetic spots are very unlikely 
in these cool atmospheres. 

At this point we cannot distinguish between a scenario 
where a high altitude cloud or storm feature is respon- 
sible over a scenario where the atmosphere is composed 
primarily of broken clouds and clearings. Further mon- 
itoring, and improved atmosphere models where hetero- 
geneities are represented in a more physically consistent 
way will be required to confirm this result. 

6.2.1. Could the variability of 2M2139 be related to binarity? 

Given 2M2139's potentially double-peaked light curve, 
and similar amplitudes of variations observed in J, H , 
and Kg bands in at least one epoch (see figure [5]) we con- 
sider here whether the observed variability could arise 
from an eclipsing binary system. Since the continuous 
shape of the light curve is inconsistent with that of a de- 
tached binary system (there is no flat out-of-transit re- 
gion), a semi-detached configuration would be required. 
In addition, subtle changes in the light curve shape from 
night to night, as well as the overall light curve asymme- 
try would require an extra variable component, such as 
a hot-spot due to mass transfer. The requirement that 
one component fill its Roche Lobe allows us to place con- 
straints on the range of primary and secondary masses 
and ages that could form a semi-detached system. For 
a given primary mass. Mi, and mass ratio q = M2/M1 
(where M2 < Mi), the Roche Lobe size of the sec o ndary 
can be computed via the formula of lEggletonl (|1983[ ) 
given by i?L/a = 0.49g-2/3/(0.6g-2/3 +ln[l + g-^/^]). 
The separation, a, is determined from Kepler's law as- 
suming a period of 15.44 hr, corresponding to a double- 
peaked light curve. The secondary Roche Lobe sizes can 
then be compared to typical BD radii from evolutionary 
models. In figure fTSl we show the secondary Roche Lobe 
size as a function of mass ratio for primary masses of 
10, 30 and 60 MjupUer- For comparison we over-plot 
typical secondary BD radii from evolutionary models, 
assuming a primary effective temperature of ^1300 K 
(note that the system age is fixed from the combination 
of primary mass and effective temperature, yielding ap- 
proximate ages of 20 Myr, 280 Myr, and 1.8 Gyr from 
lowest to highest primary mass). The system parame- 
ters for which secondary Roche Lobe filling may occur 
are limited. We find that Roche Lobe sizes are much 
larger than typical secondary radii for almost all primary 
masses and mass ratios. This scenario is possible only for 
low primary masses (Mi < 10 Mj„p) and low mass ra- 
tios {q < 0.05 — 0.1), such that a semi-detached binary 
with a 15.44 hr period would have to be a young (< 20 
Myr). double planetary- mass system. For more typical 
BD masses and ages, the binary components should re- 
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main detached, and unable to produce the observed Ught 
curve. 

Phase variations (i.e. from the day /night side con- 
trast of an irradiated companion) in non-ecHpsing sys- 
tems can also produce periodic variations. However, in- 
trinsic phase variations of an irradiated companion would 
be washed out by the brighter primary, resulting in much 
lower levels of variability than observed for 2M2139. 

A further problem with the binary scenario is that 
multi-epoch observations spanning Aug 2009 to Nov 2009 
cannot be phased together, requiring either the Aug 2009 
or Nov 2009 epoch to be offset in phase from the other 
epochs. The timing of both epochs has been verified 
in the acquisition time-stamps as well as manual log 
books, and phasing all epochs together would require 
both to be in error — an occurrence we deem to be ex- 
tremely unlikely. In order to produce a large offset in 
the expected eclipse timings by natural means, a mod- 
erat ely elliptical, precessing o rbit is required (e > 0.15, 
e.g. iRagozzine fc Wolj (|2009[) ). This could potentially 
be ac hieved if the system is undergoing Kozai oscilla- 
tions (jKozailll962f ) due to the presence of a third body 
in a highly inclined orbit to the inner pair. However, in 
order for the precession timescale to be sufficiently fast 
(e.g. months) such a third perturber would closely ap- 
proach o r exceed the criterion for d ynamical instability 
given bv lEggleton fc Kiseleval (|1995[ ). Furthermore, the 
effects of tidal friction between the inner pair would act 
to suppress Koz ai cycles over time and cir cularize the 
inner orbit (e.g. iFabrvckv fc Tremainell200'7i ). making it 
increasingly unlikely to observe a system in this configu- 
ration. 

Based on the above considerations we feel we can rea- 
sonably exclude an eclipsing binary as the source of 
2M2139's variability. 

Alternatively, like Jupiter, BDs are expected to have 
large magnetic fields. Thus it is conceivable that some 
sort of magnetic interaction between a close binary pair 
is possible without requiring physical Roche lobe filling. 
However, such a scenario does not naturally explain a 
double-peaked light curve, and it's ability to produce sig- 
nificant and highly contrasting photospheric features is 
unclear. 

6.3. The L/T Transition 

2M2139 joins SIMP0136 as the second known cool BD 
to display persistent and significant quasi-periodic vari- 
ability in the NIR, with amplitudes well above the level of 
photon noise and systematic errors. With spectral types 
of T1.5 and T2.5, and NIR colors of J - A', = 1.3 and 
J—Ks=Q.Q8, both BDs fall squarely within the L/T tran- 
sition regime (see figure [T6l) . Although a larger sample 
is needed, observations of 2M2139 and SIMP0136 sug- 
gest that large amplitude NIR variability may be more 
common at the L/T transiti on, consistent with th e cloud 
frag mentation hypot h esis o f lAckerman fc MarlevI (|2001[ ) 
and iBurgasser et al.l ()2002|) . In addition, the positions 
of SIMP0136 and 2M2139 on the color-magnitude di- 
agram — within the L/T transition regime where the 
cloud opacity decreases significantly — provides further 
rationale to prefer heterogeneous clouds as the mecha- 
nism for the observed variability. 

Overplotted on the color magnitude diagram in figure 
[16] are vectors showing the direction of observed variabil- 
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Figure 15. Secondary Roche Lobe radius as a function of mass 
ratio, q, for primary masses of 10, 30, and 60 Mj^p, assuming an 
orbital period of 15.44 lir (black lines). Secondary radii determined 
from evolutionary models for fixed q = 0.1, and assuming a primary 
effective temperature of 1300 K are overplotted (grey lines). 

ity for 2M2139 and SIMP 0136. Rather than tracing the 
approximate evolution of field BDs across the L/T tran- 
sition, the vectors point upward of the transition path, 
rising much more steeply in Mj with decreasing J — K . 
This difference refiects the fact that evolution across the 
L/T transition occurs over astronomical timescales, over 
which a significant amount of heat is lost from the upper 
atmosphere as the higher temperature T-P profile below 
a cloudy atmosphere evolves toward the lower tempera- 
tures under a clear atmosphere (e.g. iSaumon fc MarlevI 
120081 ). When we observe variability due to clouds and 
clearings at a snapshot in time we see underlying clear 
patches that are warmer than the integrated effective 
temperature of the atmosphere, and hence warmer than 
the atmosphere of a more evolved cloud-free T-dwarf of 
the same effective temperature. Therefore, the direction 
of the L/T transition path is determined by both decreas- 
ing cloud coverage and a loss of entropy from the upper 
photosphere that occurs over astronomical timescales, 
while the direction of the instantaneous variability vector 
traces only changes in cloud coverage. 

While heterogeneous clouds appear to be the most 
likely explanation for our observations, the physical 
mechanism by which clouds may fragment at the L/T 
transition remains unknown. One possibility is that as 
clouds form progressively lower in the atmosphere as the 
BD cools, eventually the entire vertical extent of the 
cloud is found within the dynamic troposphere which is 
more subject to local variations in updrafts and down- 
drafts arising from convection. These localized weather 
patterns are better able to disrupt the cloud than when 
the BD is warmer and a substantial fraction of the cloud 
opacity is found within the relatively quiescent strato- 
sphere. While weather remains a poorly understood phe- 
nomenon even among the planets in our own Solar Sys- 
tem, BDs represent a simplified case where atmospheric 
dynamics result primarily as a consequence of rapid ro- 
tation and internal heat, without the complication of ex- 
ternal forcing due to irradiation from a parent star. Thus 
future observations of weather in BD atmospheres may 
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Figure 16. Color-magnitude diagram in the MKO filter system 
for M, L, and T dwarfs with measured parallaxes (black, blue and 
red points respectively) from the compilation of S. K. Leggett. The 
approximate positions of SIMP0136 and 2M2139 are overplotted 
based on spectroscopic parallaxes as opposed to real ones. The 
J — Kg color of 2M2139 (converted to the MKO system) is based 
on the colors derived from WIRC and SpeX epochs (see table |4]l, 
as opposed to 2MASS. The approximate position of 2M2139 on 
the color magnitude diagram based on its 2MASS color (also con- 
verted to the MKO system) is shown as a green square (its absolute 
magnitude was determined relative to the WIRC/SpeX epochs, as 
the empirical relations for determining absolute magnitude have 
not changed). Solid black lines from right to left indicate the 
/sod =li2,3,4, and nc model colors respectively for log 3 =5.0. Al- 
though we have used spectral models with logg = 4.5 throughout, 
a value of logg = 5.0 is more appropriate for comparison with 
the population of field BDs (in addition, see the discussion in sec- 
tion 14.21 about choices of log g values) . Arrows depict the vector 
direction of the variability observed for 2M2139 (this paper, aver- 
aged from the 3 epo chs in figure |4]l and reported for SIMP 0136 
IIArtigau et al.ir2009l) . both magnified in amplitude by a factor of 
3. 

provide a novel opportunity to study atmospheric circu- 
lation and cloud meteorology in a higher gravity regime, 
never before probed. 

7. CONCLUSIONS 

We have reported the unique, large-amplitude quasi- 
periodic variability of the cool L/T transition BD 
2M2139. Our major findings and conclusions are sum- 
marized here. 

1. Variability We have observed 2M2139 to be highly 
variable in the NIR over 9 separate epochs span- 
ning 100 days, with an amplitude as high as 26% in 
the J-band. Measured amplitude ratios range from 
Aj/Ah = 0.84 - 0.91 and AkJAj = 0.45 - 0.83, with 
uncertainties of ±0.07-0.15. It remains unclear whether 
these variations in relative amplitude from epoch to 
epoch are intrinsic, or reflect hidden measurement un- 
certainties. Disagreement at the 2a level between epochs 
suggests the former. In addition to the short timescale 
variability, comparison with 2MASS photometry hints 
that there may be a longer-term (10 yr) variation that 



consists mainly of a brightening in the J-band. 



2. Periodicity We have measured a minimum period of 
7.721±0.005 hr for 2M2139. It is also possible that the 
light curve is double-peaked with a period of 2 x 7.72 = 
15.44 hr. In the latter case, the period would be some- 
what larger than the 2-12 hr periods typical of field BDs. 

S. Time Evolution of the Light Curve Despite a clear pe- 
riodicity, the light curve shape changes from night to 
night, and we cannot simultaneously phase our earli- 
est and latest epochs with the middle ones, indicating 
substantial evolution of cloud features on timescales of 
weeks to months, with more subtle changes occurring on 
a timescale of days. The differential rotation of a storm 
feature is a possible explanation for the longer timescale 
evolution observed. 

4- Modeling variability due to patchy clouds We have 
searched for linear combinations of ID atmosphere mod- 
els differing in temperature and condensate properties, 
that can simultaneously reproduce the observed multi- 
band variability an d the NIR spectrum of 2M 2139. Us- 
ing the models of iSaumon fc MarlevI (|2008f ) we found 
that the observations can be reproduced by a heteroge- 
neous surface wherein regions of higher condensate opac- 
ity (/sod=l or /scd=2) are cooler, and regions of lower 
condensate opacity (/scd = 3 or /scd = 4) are warmer 
(by ~ 175-425 K). The best-fitting model suggests that 
cool, thick cloud features — for instance one or more 
storms systems occupying 25% of 2M2139's visible disc 
— may be responsible. Alternatively, the data is also 
well-reproduced by a model where atmospheric hetero- 
geneities consist of warm clearings in a cloudy atmo- 
sphere, however "clear" regions must retain some degree 
of cloud opacity (i.e. /sed=4). The large temperature 
contrast required between thin and thick cloud patches 
may suggest that the thick clouds extend to higher alti- 
tudes in the atmosphere. We reiterate that these results 
should be viewed cautiously, as interpolations between 
independent ID models do not necessarily result in phys- 
ically plausible configurations. 

5. Origin of Variability The continuous, asymmetric, and 
evolving nature of the light curve strongly suggests that 
atmospheric cloud features are responsible for 2M2139's 
variability. Due to 2M2139 's (i)status as an L/T tran- 
sition BD and the expectation of patchy clouds in this 
regime, (ii)cool, neutral photosphere, and (iii)our model- 
ing of the observations we have argued that surface fea- 
tures responsible are heterogeneous clouds rather than 
magnetically induced spots. We have also considered 
that 2M2139 may be an interacting binary, but have 
ruled out this scenario. 

6. The L/T Transition Both 2M2139 and SIMP0136, 
the two known examples of BDs that display large- 
amplitude, persistent variability in the NIR, have colors 
and spectral types falling directly within the L/T tran- 
sition, suggesting that variability in this regime may be 
higher-amplitude and/or more frequent. These results 
therefore provide empirical support for the idea that the 
fragmentation of dust clouds can explain observed prop- 
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erties of the L/T transition. 

Looking forward, long term monitoring of this target 
both photometricaUy and spectroscopically, and over a 
broader wavelength regime, should reveal the true nature 
of its variability. In addition, a parallax measurement 
and high resolution speetrum will be essential in order 
to better constrain 2M2139 's absolute magnitude and 
physical properties. 
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APPENDIX 

A. RELATIVE FLUX CALIBRATION OF SPEX PRISM LIBRARY SPECTRA AND SYNTHETIC 2MASS COLORS 

The relative flux calibration between J, H , and Kg bands for M, L, and T dwarf sources in the SpeX Prism Library 
was investigated by determining synthetic 2MASS colors from the SpeX spectra and comparing them to reported 
values from the 2MASS catalog. S ynthetic 2MASS col ors were found using the relative spectral response curves and 
zero-magnitude fluxes provided bv I Cohen et al.l ()2003[ ). The error estimates, aspeX are the sum of a random error 
component determined from the measurement uncertainties associated with individual wavelength bins (typically 
small), and an additional 0.07 mag corresponding to the average difference in synthetic J — Kg colors computed for 
(presumed non-variable) targets which have been observed at two epochs within the SpeX Prism Libraries. 

Differences between 2MASS colors and SpeX synthetic colors were determined for all M, L and T dwarf sources 
from the SpeX Prism Library that also had moderate to high quality detections in 2MASS (2MASS catalog QFLG of 
at least 'C'). For each BD meeting these requirements, the difference between the 2MASS and synthetic SpeX colors, 
A( J — Ks), was computed and assigned an uncertainty of 

where aj and ax are the photometric errors taken directly from the 2MASS catalog. The resultant A{J—Ks)/aj-K^ 
for all BDs considered is plotted in figure[T71 and shows that 2M2139 is a 4cr outlier. If we consider only L and T-dwarfs 
the match between 2MASS and SpeX colors is surprisingly good; a gaussian fit to a histogram of A(J — Ks)/aj^K^, 
also shown in figure [T71 has width of cr ^1 and is roughly centered about zero (//=-0.14). While differences in J — Kg 
color are roughly symmetric about zero for L and T dwarfs, there appears to be a small systematic offset for M dwarfs, 
the origin of which is not clear. 

Therefore, uncertainties in synthetic J — Kg colors for SpeX Prism Library L and T dwarf spectra are well estimated 
by propagating the standard errors reported for each wavelength bin, and then adding to this 0.07 mag. 

B. MODELS FOR 2M2139 S VARIABILITY AT ADDITIONAL EPOCHS 

In section [Ql we presented simultaneous model fits to both the NIR spectrum and multi-color variability of 2M2139. 
The model fits presented were constrained using Ah I Aj and Ak^ / Aj amplitude ratios measured from the simultaneous 
JHKg light curves obtained on 30 Sep 2009 (figure [T^ . Here we provide equivalent model fits constrained using 
amplitude ratios measured from the other two epochs for which we obtained simultaneous JHKg light curves: 26 Sep 
2009 (figure nil), and 01 Oct 2009 (figure [HI). 
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Figure 17. Comparison between 2MASS photometry and synthetic 2MASS colors derived from SpeX prism spectra for all L and T dwarfs 
in the SpeX prism library. Sources flagged in the prism library as having a poor SNR, or those with poor quality J and/or Ks 2MASS 
photometry (Point Source Catalog quality flag not equal to 'A','B',or 'C'), were omitted. 
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Figure 18. Same as figure [T2l but with amplitude ratios constrained using data from the 26 Sep 2009 epoch. 
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Figure 19. Same as figure [T2l but with amplitude ratios constrained using data from the 01 Oct 2009 epoch. 
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